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Abstract 

We construct a matrix model equivalent (exactly, not asymptotically), to the random 
plane partition model, with almost arbitrary boundary conditions. Equivalently, it is 
also a random matrix model for a TASEP-like process with arbitrary boundary 
conditions. Using the known solution of matrix models, this method allows to find 
the large size asymptotic expansion of plane partitions, to ALL orders. It also allows 
to describe several universal regimes. On the algebraic geometry point of view, this 
gives the Gromov-Witten invariants of with branes, i.e. the topological vertex, in 
terms of the symplectic invariants of the mirror's spectral curve. 



1 Introduction 

The statistical physics problem of counting plane partition configurations of some do- 
main, as well as its various equivalent formulations, has become a very active and fas- 
cinating area of mathematical physics in the past years, culminating with Okounkov's 
renowned work. Beyond a beautiful combinatorics problem, it has also many indirect 
applications, like a tiling problem similar to a discrete version of TASEP, i.e. the sim- 
plest model of out of equilibrium statistical physics, and algebraic geometry, as it plays 
a key role in the computation of Gromov-Witten invariants of some toric Calabi-Yau 
3-folds, through the topological vertex method [3]. 

The works of Okounkov, Kenyon and Sheffield [IB], have brought immense progress, 
in the understanding of large size asymptotics behaviors of plane partitions. It was 
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observed, that in many universal regimes, the statistical properties of large plane par- 
titions, is very similar to that of matrix models, and many works have taken advantage 
of that similarity. 

Here, in this article, we show that there is not only a "similarity" between plane 
partitions countings and matrix models, in fact we show that plane partitions IS a 
matrix model, even for finite size. As a consequence, we may use all the machinery 
developped for solving matrix models, and we are able to compute all orders corrections 
to the large size asymptotics. 

Our matrix model, is a multimatrix model, with non-polynomial potentials. It may 
look very complicated at first sight, and its spectral curve may look rather compli- 
cated too. However, the solution of matrix models is expressed in terms of symplectic 
invariants, and up to a symplectic transformation (which does not change the sym- 
plectic invariants), our complicated matrix model's spectral curve, is equivalent to the 
Harnack curve of Kenyon-Okounkov-Sheffield [18] . 

Moreover, our formulation allows to use the full toolbox of matrix models technol- 
ogy. For instance the method of orthogonal polynomials gives determinantal formulae 
for correlation functions, the integrable structure, Riemann-Hilbert problem, and much 
more. And the loop equations method allows to compute the large size expansion order 
by order [SU [53] . 

1.1 Main results of this article 

Our main result is the theorem l3.lt 

Theorem I3.lt The tiling model-plane partitions-tilings generating function, can be 
written as a matrix integral, (a more precise result is written in theorem \3. This 
identification is exact, it is not asymptotic. 

An immediate consequence is obviously: 
Corollary: All asymptotic limits of tiling model-plane partitions-tilings, are random 
matrix limit laws. 

(but it remains to classify all possible random matrix limit laws). 

Also, since our matrix model is a chain of matrices, classical results of matrix models 
apply: 

Corollary: The matrix model is integrable, the generating function is a Tau-function, 
and for instance correlation functions are given by determinantal formulae of Janossi 
densities type f31f . 

We would like also to emphasize that our "Tiling ^ matrix model" identification 
works for very general cases, with almost any possible boundary conditions, we can 
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also give weights to points of the domain to be tiled, some points can be forbidden 
(defects), or obliged, or just have an arbitrary weight. 

The second part of our paper, starting at section |U consists in "solving" the matrix 
model. We do it explicitly only for some not too complicated boundary conditions. We 
recover the Harnack curves of Kenyon-Okounkov-Shefiield. 

We discuss many examples in sections [7] and after, and in particular we apply our 
method to the enumeration of TSSCPP's in section [HI 
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2 Some statistical physics models 
2.1 Plane partitions 

Consider 3 integers Nx, N^j^, N^, and 3 partitions A, /i, z/, for example: 

A = /x = t z/ = fc (2-1) 

A plane partition vr with boundaries A, /i, u, and of size A^a, A"^, A^^, is a piling of cubic 
boxes in the corner of a room, with boundary conditions given by A, /i, u, for example: 




A^A is the height of the plane-partition (height of the cubic boxes piling), A^^ (resp. N^) 
is the extension towards left (resp. right), so that beyond A"^ (resp. N^,), the section 
is frozen to fi (resp. z/). 

The partition function we would like to compute is: 

TT, dlT=(X,fl,u) 

where |7r| is the number of boxes, i.e. the volume, it is called the "weight" of tt. 

This partition function is the so-called "topological vertex" in topological string 
theories [3], it is the building block to compute Gromov-Witten invariants of all toric 
Calabi-Yau 3-folds [3[Tl[3l[52l[56l[571[6ll[62l|67l[55]. 

From the combinatorics point of view, it is the generating function for counting 
plane partitions with given boundaries and weighted by their volume. From the sta- 
tistical physics point of view, it can be viewed as a model for a growing 3-dimensional 
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crystal in the corner of a room. All those topics have remained important research 
areas in physics and mathematics, and it would be difficult to summarize all what has 
been done. Let us mention that Kasteleyn |15] found an explicit expression for the 
partition function of a domino-tiling, which can be rephrased as plane partition, and 
since then, the subject has been studied a lot, see for example [IHl SZl | 



2.1.1 Remark: semi-standard tableaux 

If we slice our plane partition tt at all integer times (time = horizontal coordinate) 



-N, 



.N^, at each time the slice is a 2-dimensional partition A(t). 



-2 



-1 



-5-3 



1^ 



_Q 



-2 



-1 







V 



>6 



_L 



Figure 1: A plane partition is equivalent to the data of two semi-standard Young 
tableaux of the same shape A(0). Each of these two tableaux, is the superposition of 
growing (or decreasing) partitions A(t). 



It is then clear that partitions A(t) are growing from t = — iV^ to 0, and then 
decreasing from t = to N,^: 



Wt<0, A(t)^A(t + l) 



Vt>0, A(t)^A(t-l) 



(2-3) 



where -< is the partial ordering of partitions (A -< /i means that A can be obtained from 
IJ, by removing boxes). Since we have /i -< A(— A''^) ^ . . . X{t — 1) ^ X(t) -<...-< A(l) -< 
A(0), we may draw all A(t) with t < inside the Ferrer diagram of A(0), and we write 
in each box, the time t at which the box appears for the first time. We do the same 
for t > 0, and we have two semi-standard tableaux with the same shape A(0). A semi- 
standard tableau is a Ferrer diagram with integer entries decreasing along columns, 
and strictly decreasing along rows. See figUl 
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When A = /i = z/ = 0, and N\ = Nf^ = = oo, the statistics of the partition 
A(0) is the sum over all pairs of semi-standard tableaux of shape A(0), i.e. it is the 
Plancherel measure [Ml ESl [69] : 



'^('^) = ( I , [, I = (#(semi — standard tableaux of shape A))^ (2-4) 

V 1^1' / (1^10 

We shall study this limit in section 17.41 
2.2 Jumping non- intersecting particles 

T.A.S.E.P. means "totally asymmetric exclusion process" [51], it is the simplest model 
of statistical physics out of equilibrium, it has focused considerable amounts of works 
[TSl IT9l [63l l35l [36l [38] . and it is still intensively studied. It is a model of self avoiding 
particles which can either stay at their place, or jump 1 step forward, provided that 
the next space is unoccupied. In the dynamics we shall be considering here, time is 
discrete, and at each unit of time, several particles can jump. 

It is well known that plane partitions can be rephrased in terms of self-avoiding 
jumping particles model which is a kind of discrete T.A.S.E.P. [Ml [65| [69l [68| [^ [6| [5| [2] . 
let us re-explain it here. 

Let us draw in red the A^a non-intersecting lines going through tiles and 

The tiles correspond to upright lines with slope +1/2, and the tiles ^ correspond 
to downright lines with slope —1/2. The intersection of those red lines with integer 
time hues t = —N^, . . . .N^ are interpreted as positions of some particles hiit) {h = 
is at the top). See figure [2l 

Therefore, a plane partition vr, can also be described as the data of 1 + A"^ + A"^ 
sets of A'a variables: 

hiit) , te{-N^,...,N,} , t = l,...,Nx (2-5) 

such that for every integer t G {— A"^, . . . , N^,}, the hi{t) — are non-negative ordered 
integers: 

hi{t)-^\t\eN , h{t)> h2{t)> h{t)> ...> hN^{t)>^\t\+rt{X) (2-6) 

where rt{X) is the profile function of the partition A. 
At time t = — A^ and time t = we have: 

hi{-N^) = + Nx + ]^N^ , hi{N,) = Vi-i + Nx + ^N, (2-7) 
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Right 



1^ 




Figure 2: The level lines, go only through upright or downright tiles, they form A^a 
non- intersecting lines whose slopes are piecewise 



and at each time t, we have a partition \(t) = {\i(t) > \2{t) > • • • > '^Nxi't))- 

\i{t) = h,{t) + z - Nx - ^\t\ 



(2- 



we thus have A(— A*"^) = fi and \{+Ni, 
Moreover we have: 



hi{t) - hi{t+l) = ±1/2 



(2-9) 



This is what we call here a discrete T.A.S.E.P. process: 

• there are N particles at positons hi{t), and at each unit of time, they jump by 
±|, and they can never occupy the same position (in the usual formulation of TASEP, 
particles jump by or 1, and here we have tilted the picture so that they jump by 
±1/2, which is clearly the same thing up to hi{t) hi{t) + 1/2). 

2.2.1 Summary plane partition and jumping particles 



A plane partition tt with boundaries A,/i, i/, and of size Nx, N^, N^, is equivalent to 

it: 

hi{t) -hi{t+ 1) = ±1/2 



hi{t), i = 1, . . . , Nx, t = —Nfj_, . . . , A^;^, such that: 
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• hi{t) > h2{t) > hit) > ...> hN,{t) > i|t| 

• hi{N,) = Ui - i + Nx + ^N, 



The total numer of boxes in the partition is: 



\TC\ 



(2-10) 



(2-11) 



2.2.2 Height function and density 

There is a relationship between the height function of the pile of cubes, and the density 
of /li's. Define the density at time t as the Dirac-comb distribution: 



1 

^ 1=1 



(2-12) 



The profile of the partition at time t is recovered from the integral of p as follows. 
Define the integrated density: 

I{h, t) = Nx + Nx / p{h', t) dh\ (2-13) 



I{h,t) computes the index i such that h{x,t) = hi. Then, define the function 



X{h,t) = h - -\t\+ I{h,t) - Nx. 



(2-14) 
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The plot of X{h,t) against I{h,t) is the shape of the partition is X{t) at time t. See 
figJl 




Figure 4: The density p{h, t) = ^ S{h — hi{t)) encodes the profile of the partition 

A(t). 



The surface of the pile of cubes in 
at all times, i.e. it is given by: 



is recovered from plotting the partition A(t) 



a;i = A + -t) = I -Nx + h + 
X2 = X + l{\t\ +t)=I-Nx + h- 

X3 = I 



(2-15) 



2.3 Lozenge tilings 

Another representation of plane partitions and the self-avoiding particle model, is with 
lozenge tilings of the rhombus lattice. 

The rhombus lattice is a tiling of the plane, with lozenges ^C>, whose centers are 
at the positions {h, t) such that /i — | G Z and t G Z, see fig El 

A plane partition or a self-avoiding particle model configuration can be represented 
as = A^A oriented self-avoiding walks in some domain T> of the rhombus tiling of the 
plane. See figure [61 At a given time, each walker moves to the right t —>■ t -\- 1, and 
either np h h — ^ or down h ^ h + ^. 

The centers of occupied lozenges are at positions {hi{t),t), i = 1, . . . , N. 

The domain V is the domain which contains all possible paths. P is a subdomain 
of the hexagon defined by the following 6 inequalities: 



< 



t 



< U 
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Figure 5: The rhombus lattice is a tihng of the plane, with lozenges whose centers are 
at the positions {h, t) such that /i — | G Z and t G Z. Notice that we orient h from top 
to bottom. 



t ^min ^ 7 ^ 1 I !_ \ I ^ ^min 



hN {tinin) 7: ^ h < /ll(tmin) + 

Z z, 

'^max i ^ 1 ^ 1 { ± \ , ^max ^ 



/i7v(W) 7^ < h <hi{t 

(2 - 16) 

Indeed, since at each time, and particularly at t = tmin and t = tmax, we have /iat < 
hi < hi, no path starting at some /li(tmin) at t = tmin and ending at some hj{t.^a_^ at 
i = ^max can go out of this hexagon. 

Moreover, the domain D can be chosen such that some positions are forbidden, in 
fact we allow "D to be any arbitrary subset of the maximal hexagon. 

For example for plane partitions, the positions corresponding to the boundaries 
A,/i, 1/ are forbidden, i.e. V has holes corresponding to the 3 partitions A,/i, z/ at the 
boundaries: 

• At t = tinin = we remove from T>, all lozenges which are not of the form 
/li(yu) for some i. 

• At t = tjnax = A^j,, we remove from P, all lozenges which are not of the form hi{v) 
for some i. 

• At every time t we impose /i > y + rt(A). 
See fig El 

2.3.1 Defects 

It is interresting to generalize our model to include walks on more general domains, 
not only limitted by 3 boundary partitions. In particular, we may allow defects and 
holes at almost any place. 
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Figure 6: A plane partition configuration, can be represented as self-avoiding walks 
in some domain of the rhombus tiling of the plane. The white lozenges are forbidden, 
they are related to the 3 boundary partitions A, /i, v. 



Consider a connected compact domain V in the rhombus lattice, but not necessarily 
simply connected. 
Definitions: 

The maximal domain Max(I') of T> is the intersection of tmin < t < tmax with 
^Af(imax) + < h < /ii(tmax) ^ (see the dashed region in figure^ and figure 

The defect ofD, isT) = Max('D)\T'. The shadow ofD, is the domain V C Max(T') 
which is unaccessihle to particles moving in T>, with slopes We have V (ZT>. 
Those notions are illustrated in fig [HI 

Definition: 

The minimal defect Vq ofV, is the smallest suhdomain ofV, such that: 

Vo = V (2-17) 

2.3.2 Domain at time t 

Let us call the slice of the domain T> at time t, i.e. the set of allowed positions 

of particles at time t: 

V{t) = {h/ h-^eZ, {h,t) eV} (2-18) 
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I ' t ^ 

min max 



Figure 7: The maximal domain of T> is the region comprised in the dashed hne. 



Let us call 'D{t) the slice of V at time t, i.e. the position of holes: 

V{t) = {h/ h-^eZ, {h,t) eV} (2-19) 
Let us assume that the number of particles N, is such that: 

Vt ^miri; • • • ; tmax , #^^(t) > N, (2-20) 

in other words such that it is indeed possible for avoiding particles to move in V. 

Remark 2.1 Most often we will choose N = #I'(tmax)) i-e. the position of all particles at 
time tmax are fixed. Very often we will also choose the domain such that N = #'D(tmax) = 
#D(tmin), i.e. the position of all particles are also fixed at time tmin- 

But we emphasize that the method we present here works also without those assumptions. 

2.3.3 Filling fractions 

In general, the domain at time t, is a union of intervals: 

mt 

V{t) = U [ai{t)Mt)]- (2-21) 
14 




Figure 8: Consider a given domain V. Its maximal domain Max(D) is represented with 
dashed hne on the left figure (it is obtained by lines of slopes ±| starting from the 
extremities of the domain at tmax)- In the middle figure we have represented the 
defect of that is "D = Max(P) \ V (the red region). In the right figure we have 

represented the shadow of V, that is the domain V C Max(T') which is unaccessible 
to particles moving in V, with slopes ib| (the brown+red region), it contains V. This 
means that particles moving in T) can never enter the brown tiles. 



If the intervals are disconnected, no particle can jump from one interval to another, and 
thus the number of particles moving in an interval is constant in time, until intervals 
join. We may thus fix the number of particles in each interval. Call it ni{t) = filling 
fraction of the interval [a^, 6j] at time t. We have: 

mt 

n,{t) = N. (2-22) 

i=l 

In general, if the domain has k holes, there are k independent filling fractions. 

2.4 Partition function 

The partition function which we wish to compute is: 

+ _i 

t-max 2 

Zj,{V;q,a,(3)= J] gE,./^.W JJ a(t')*(^)(*') /?(t')*^^^'*'' 

(2-23) 

where at each time t we have hi{t) > h2(t) > ... > h]^{t), and where /^^(tmm) and 
hi{tma.x) have fixed values given by two partitions /i and u (encoded by the boundaries 
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Figure 9: The left figure represents tfie domain T> and its defect T> in red. The right 
figure represents (in red) the minimal defect Vq of P. The middle figure shows that 
both V and have the same shadow (in brown+red). 



of V): 

hiitrain) = hi{fi) , hi{t^^^) = hi{v) (2-24) 

We count the configuration with a weight g''^' oc g^*-' and with weights a(t + 1) 
per number of upward jumps between time t and t + 1, and with weights (3{t + |) 
per number of downward jumps between time t and t + 1. 

It can be rewritten: 

ax 

Z^{V-q,a,(3) = n n 

hl{t)>...>hN{t), h^{t)eV{t) i = l t=tynin 

n nH^')M^^(^'+2)-^^(^'-2^-2) 

t'=tmin + l/2 i = l ^ ^ 

(2 - 25) 

where here, 5 is the Kroenecker's (5-function. 
2.5 Applications 

• In topological string theory, Gromov Witten invariants of toric Calabi-Yau 3-folds are 
computed with the topological vertex, which is the following sum of plane partitions 
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[21 EH]: 

Zoo(g,l,l;A,/i,z/)= Yl (2-26) 

7r,97r={A,^,!^) 

• When X = = u = (/}, and N\ = = iVy = cxd, and a = P = 1, this sum is known, 
it is the Mac-Mahon formula: 

oo 

Zoo(g,l,l;0,0,0) = J] g'"' = JJll = 1 + 9 + 3g' + + 13g^ + .. . (2-27) 

TT A;=l 



• The Razumov-Stroganov conjecture [U [22] has put forward a problem of combi- 
natorics of totally symmetric sel-complementary plane partitions (TSSCPP), and the 
claim is about the relationship with the combinatorics of alternating sign matrices 
(ASM). There is a 6- vertex matrix model formulation for ASM [72], and it would be 
interesting to also have a matrix model for TSSCPP, that's what we address in section 
|8] below. 



3 Matrix model 



We are going to represent our self avoiding particle process partition function eq. (12-231) 
as a multi-matrix integral. 

Let us sketch the idea of the next subsection: we shall introduce tmax — ^min + 1 
normal matrices Mt of size N x N, for all integer times t between tmin and tmax, whose 
eigenvalues are the hi{t). Moreover, we shall Fourrier-transform the 5— functions which 
enforce hi(t + 1) = hi{t) ± 

oi{t + ^) 5{K{t + 1) - hiit) - ^) + Pit + 1) 5{K{t + 1) - Kit) + 1) 

+ 00 



(3-i; 



oo 

+ 00 



a{t + 



^\ 2i-Kn{h,(t+i)-h,{t)~\) 



+ /?(t+ -)e' 



2iTTri{hi{t+l)-hi{t)+\) 



c(t + l)e-^-'-'+/3(t+l)e^--» 



i.e. we shall introduce some Lagrange multipliers rj(t + 1) at all half-integer times, and 

whose eigenvalues are the Lagrange 



we will introduce a N x N hermitian matrix R,, i 

^^2 



multipliers ri(t + | 



2/' 



1, . . . , A^, which implement the 5— functions for the jumps 



between time t and t + 1. See figJTOl 

More generally, if we wanted to allow jumps of several steps hi{t + l] 
{si, . . . , Sfe}, we would take the Fourrier transform of aje^*'^*^^*. 



hi{t) + s, s e 
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The non-intersecting condition for paths can be reahzed as a determinant, hke 
Gessel-Viennot formula [37], which allows to rewrite 

Jj-g2i^r,/., ^ det(e2^""'''0 (3-2) 

i 

and we recognize that this expression is the Itzykson-Zuber-Harish-Chandra formula 
[391 iU: 

det(e2^""^'^0 = A(i?)A(/i) / dU (3-3) 

Ju{N) 

This is the key to obtain a matrix integral, it introduces angular degrees of freedom in 
addition to eigenvalues: 

Mt = Uh{t)U^ , /i = diag(/ii(t),/i2(t),...,/i7v(t)). (3-4) 




3/2 5/2 7/2 t+1/2 



Figure 10: We introduce a chain of matrices. The eigenvalues of the random matrices 
Mt with t integer, are the random hi{t). And the Rt' with t' half integer are Lagrange 
multipliers which enforce the relations hi{t + 1) — hi{t) = 



And we shall introduce some potential Vt for the matrices Mj, to ensure that their 
eigenvalues are in X'(t), and Vt^^^^^ and Vt^^^^ enforce the initial values at time t = tmin 

or t ^max- 

We shall find that the sum over /ii(t)'s can be rewritten as a "chain of matrices" 
matrix model. 

So, let us describe the model now. 
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3.1 The multi-matrix model 



Let us consider the following multi-matrix integral: 



tmax — 1 /. " — ""lax j 

J D 



Z = W dMt W dR 

^—^min t-ti-nin+2 

imax "-m^ix j ''max j 

Tr Vt(Mt) ^TrMt TT Tr f/,, (ij,, ) TT „ Tr i?^, (M^,_^ i -Af^,_ i ) 



J-j- g-Tryt(Mt) ^TrMi J-j- g- Tr f/,, (ij,, ) J-j- 



"2 ^ ~2 

*=tmin t'=tra\a+\ t'=imin+5 

(3-5) 

Each integral over Mt with t = tmiru • • • ? ^max — 1 is over the set Hj^ of hermitian 
matrices of size N , and each integral over Rti with t' = tmin + I, • • • , ^max ~ | is over the 
set iHfyf of anti- hermitian matrices of size A^. And there is no integration over Mt^^^, 
which is a fixed external field, which we choose equal to: 

Mw= diag(/ii(z/),...,/i;vM) (3-6) 

The potentials Vt or Uf are defined as follows: 

• for t' half-integer: 

e-^*'M = «(t')e-i+/3(Oe5 (3-7) 

this potential for Rf is the Fourrier transform of the jumps, see eq. (13-11) . we have just 
rescaled r by 2i7r. 

• and for t integer, t^ain < t < tmax, we choose the potential Vt{x) such that: 



t „ _ ( iifxev{t)/v 

•2 



Vx G (- + Z) n Max(I)) , e-^*(^) = <^ o if a; G Vo{t) (3- 



arbitrary otherwise 

e~^* can be more or less interpreted as the characteristic function of the domain P(t), 
or more precisely, the complementary domain of the defect T>{t). 

However, this definition does not define a unique potential Vt, and many potentials 
Vt may have the property eq. (13-81) . In particular, we see that the value of Vt{x) in the 
shadow of V, more precisely in V/V^, is undetermined. We show below some rather 
canonical examples of Vt satisfying those constraints. 

We will see below that the value of the partition function Z does not depend on 
the choice of Vf. 

• and for t = t^in, we choose the potential Vt^.^{x) such that: 

, ( ^0 if X = hi{fi) 

Vx G (^ + Z) n Max(P) , e-^*-i"(^) <^ = if x ^ hi{ix) (3-9) 

I arbitrary otherwise 
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In other words, we do not even require that e ^'min = 1 in the domain, we only require 
that it is 7^ 0. 

Examples of Vt: 

Since our domain is compact, we have to find a potential Vt with prescribed values 
at a finite number of points, and a possibility is to choose e"^* to be the Lagrange 
interpolating polynomial going through the prescribed values. 

For example, if our domain at time t is contained in the interval 

V{t) C[kt-,kt+], (3-10) 
we may choose the potential Vt as the Lagrange interpolating polynomial: 

e-v-.W = l_ Y: . \pr^^\. ' Pk-,H{h)=ll(h-j). (3-11) 

Another possibility is to take the limit fc+ — *• cxo: 

- y (3-12) 



i\ [x — i) r(— x) 



or we may also choose: 



e 



-^*(^) = y (3-13) 



Depending on the type of applications we are interested with, it is sometimes more 
convenient to work with a potential of type eq. (13-111) or a potential of type eq. (13-121) . 
or also their g-deformations, or sometimes other potentials having property eq. 



In general, we see that —Vt{x) must have logarithmic singularities on Po(^), and 
thus we may write: 

Vl{x) = -Y.^- +f'i^)+9'{x-l) (3-14) 
X — t 2 

where g{x) can be any arbitrary entire function which vanishes on x G Z, and where 
f{x) is an analytical function such that if i G T>{t) we have Vt{x) = 0. Since fixing the 
values of / and g does not fix the values of /' and g' at those points, we see that what 
characterizes VJ' is that it has simple poles with residue —1 in 'Do{t), plus an almost 
arbitrary analytical function. 
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3.1.1 Diagonalization 

Let us diagonalize all matrices in eq. (13-51) . One needs to know that any normal matrix 
(and in particular hermitian matrices) can be diagonalized by a unitary transformation: 

Mt=UtXtUl , WiGf/(iV) ,X, = diag(Xi(t),...,X^(t)) (3-15) 

and the matrix measures are: 

dMt = ^ dUt dXt ^{Xtf (3-16) 
where dlAt is the Haar measure on U{N) and A(Xf) is the Vandermonde determinant: 

A{Xt) = lliX,it)-X,{t)) (3-17) 

i>j 

Similarly for the matrices Rt': 

Rt> = Ut' Ye Ul , Ut, G U{N) , Ye = diag(ri(t'), • • • , ¥^{1')) (3-18) 
Therefore we may rewrite eq. (13-51) as: 

J^jY1^2(tmax-tmin) + l ^ 

._. Jim , 1 

, 1 

tmax-1 'max 5 



max 



i=imiii t'=t ■ +- 



2 



Yl I{X,_,.,Ye)I{-X,^.,Ye 

f'=f ■ 

'' ''min~ 2 

(3-19) 

where /(X, Y) is the Itzykson-Zuber integral 



I{X,Y)= / ciWe^^^"^"' (3-20) 

IU{N) 



It is well known that [391 

T(X Y\ = 

A(X)A(y) 



nx,y)= 'ri^:!!:! (3-21) 



and therefore [53], we get: 

J^jY1^2(tmax-tmin) + l ^ 
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— 'rnin =T. :„ -I 



ff=f ■ A-L 



tu 



TrC/t/(i;/) 



1_ 

J]' det(e^'(*'+5)^.(*')) det(e-^«(*'-5)^^(*')) 
(3 - 22) 

3.2 Relation between the matrix model and the self-avoiding 
particles model 

3.2.1 Integrals over Lagrange multipliers 

Then let us perform the integral over Yj[t'), we have: 



dYt, det(e^^(*'+^)^^(*')) det(e-^'(*'-^)^^(*')) e''^^''^^^*')) 



(3 - 23) 



+/5(t')5(X.,(.)(t' + i) - X,^,(,)(t' - ^) + ^) 
iV! det [a(t')5(X,(t' + 1) - X,(t' - ^) - ^) 



This term implies that there must exist some permutation such that 

^.(i' + ^)=^..w(^'-^)±^ (3-24) 

with respective probabilities (3{t') , a{t') . And in particular, since this is true for t = 
^max5 we have Wt: 

-\ezr\ [/^;v(Wx) - /^i(tmax) + ^^^^] (3-25) 

In particular that implies that Xi{t) G (| + Z) n Max(r'), and thus e"^'(^'W = 1 if 
Xiit) e "D and e~^'*^^'*^*^) = if Xj(t) e T^oif). In particular, the matrix integral is 
indeed independent of the choice of V^, provided that it satisfies eq. 
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Then, since the quantities we are summing are symmetric, we can assume, up to a 
multiphcation by A^!, that all Xi{t) are ordered: 

Xi{t) > X2it) > ...> X^it) (3-26) 

and this yields a factor (iV! )*'"=''' which we shall discard because we consider the 
partition function up to global trivial constants. 

In other words, the result of the integral is a sum over Xi{t) = hi{t), where hi{t) — | 
are ordered integers: 
1 



(t) 



^ ' /l,(t),i = l,...,Af,t=tmin + l,-,tmax-l 

n \^^^^K{t'+\)-h,{t'-\)-\^ i^{^)Ki{t'^\)-K{t'- 



-^+- 
i,t' 

(3 - 27) 

where the sum over hi{t), is such that: 

hi(t) - ^eZ , and hi{t) > h2{t) > ...> hN{t) (3-28) 

and 

/i,(W) = K{v) (3-29) 

Then, notice that e-^*('*'W) = 1 if Kit) e V and e'^'^'^'W) = if hi{t) e Voit). In 
other words, the sum is only over hi{tys such that 

hi{t) -\i Vo{t) (3-30) 

In other words, we have a self-avoiding particle process on the rhombus lattice, which 
avoids a prescribed domain V, i.e. we recover our self-avoiding particles partition 
function: 

Theorem 3.1 The self- avoiding particles model partition function Zi\f{'D,q, a, (3) in a 
domain V, is proportional to the matrix integral Z: 

Z = ^e-^«^v.(^.(.)) Z^(P,g,a,/3) (3-31) 

where the constant Cn^t depends only on N and T , and nothing else. It contains the 
normalization factors, such as the volumes of unitary groups, and the powers of N\. 

This theorem implies immediately, as a tautology, that whatever limit we consider 
(for instance large size limit, q ^ 1 limit, bulk regime, behavior near edges, . . .), the 
asymptotic statistical properties are always matrix models limit laws ! 

This explains why one finds sine kernel laws in the bulk, Tracy- Widom laws near 
some boundaries, Pearcey laws, and many more... 
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3.3 Determinantal formulae 



Theorem 13.11 allows to apply to the self-avoiding particles model and plane partitions, 
all the technology developed for matrix models, in particular the methods of orthogonal 
polynomials . 

Consider a "chain of matrices" integral: 

Z = [ f\dM, e-«T'-E-i^'(^^')+^«'=^^'^^+il (3-32) 

where Vi are some potentials, where H^iCi) is the set of N x N normal matrices having 
their eigenvalues on contour Cj, and where Mp+i is not integrated upon. 

The Eynard-Mehta theorem [31] shows that correlation functions of densisities of 
eigenvalues pi{x) = Tt 5{x — Mj), are determinants. Namely, there exist some kernels 
Hij{x,x') such that: 

(pkA^i) . ■■pkA^n)) = det , (3-33) 

where the kernels Hi j are Christoffel-Darboux kernels for some families of biorthogonal 
polynomials. 

From the point of view of self-avoiding particles model and plane partitions, this 
should allow to recover many determinantal formulae in the TASEP literature, see for 
instance [501191159]. 

Another consequence of the "orthogonal polynomials method" [M], is that matrix 
integrals of type 13-32] are Tau-functions for the Integrable Toda hierarchy pLj. This 
should allow to recover many differential equations in the TASEP literature. 

4 Matrix model's topological expansion 

The good thing about theorem 13. 11 is that the general expansion of matrix integrals of 
the chain of matrices type is known to all orders. 

4.1 Generalities about the expansion of matrix integrals 

See appendix [C] for a more detailed description. 
Consider a "chain of matrices" integral: 

Z = [ f\dM, e-'3TrELiV.(M,)+E,c.M.M.+i] (4_^) 
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where Vt are some potentials, where H^iCi) is the set of N x N normal matrices having 
their eigenvalues on contour Ci, and where Mp+i is not integrated upon. 

In some good cases (depending on the choice of potentials Vi and paths Ci), such 
an integral has a large Q expansion of the form: 



Such an expansion does not always exist. It exists only if the paths Ci which support the 
eigenvalues, are "steepest descent paths" for the potentials Vi (see e.g. Find- 
ing the steepest descent paths associated to given potentials is an extremely difficult 
problem. 

Fortunately, many applications of random matrices regard combinatorics, i.e. they 
are formal series in some formal parameter, and very often, the corresponding so- 
called "formal matrix integrals" do have a large Q expansion almost by definitioE§|, 
and eq. (14-21) holds order by order in the formal parameter (a formal parameter which 
is not necessarily Q). Here, we are considering applications to statistical physics, our 
partition functions are formal series, and we shall assume that such an expansion exists 
(order by order in a suitable formal parameter). 

The problem is then to compute the coefficients J-'g. 

The answer was found in [33], by using loop equations (i.e. Schwinger-Dyson equa- 
tions in the context of matrix models), and which just correspond to integrations by 
parts. 

The solution proceeds in two steps (which we explain below): 

1) Compute the "spectral curve" S of the matrix model. The spectral curve S = 
[x, y) is a pair of two analytical functions x{z),y{z) of a variable z living on a Riemann 
surface. The spectral curve is obtained from the "classical limit" of the integrable 
system whose tau-function is the matrix integral. Roughly speaking, if we eliminate 
z, the function y{x) is more or less the equilibrium density of eigenvalues of the first 
matrix of the chain. We explain in appendix O how to find the spectral curve of a 
general chain of matrices. We emphasize that associating a spectral curve S = {x, y) 
to a given matrix model, is something already done in the matrix models literature. 

2) Then compute the symplectic invariants Fg{S) of that spectral curve (symplectic 
invariants of an arbitrary spectral curve S were first introduced in P7], they are rather 
easy to compute, and we recall their definition in appendix [B]), and the main result of 

^For formal matrix integrals, the integration paths Ci for eigenvalues, are most often not known 
explicitely, they can be determined so that a large Q power series expansion does exist. 



oo 




(4-2) 



9=0 
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[33] is that Tg = Fg{S), i.e.: 

oo 

In Z = Y,Q'~''F'g{S) (4-3) 

9=0 

4.2 Spectral curve of the self-avoiding particles matrix model 

We recall in appendix O the main results of |33j, i.e. how to compute the spectral curve 
of an arbitrary chain of matrices. Here, in this section, we merely apply the general 
recipe of [33j (see appendix [C]) to our matrix integral eq. (13-51) . and we give a "ready 
to use recipe" . 

Recipe for finding the spectral curve of the matrix integral eq. (I3-5P : 

• Find 2(tinax~^min + 1) analytical functions of a variable z {z belongs to a Riemann 
surface C). There is one such analytical function for each matrix of the chain, plus one 
additional function at the end of the chain. Let us call them: 

('^' ^) ' ^ ^miri) • • • ) ^max ; ^(^) t ) i t '^min 2 ' ' ' ' ' ^max 2 ' (^~^) 

Those functions are completely determined by the following constraints: 

1. Those functions must obey the following system of equations \f z: 
X{z, t' + \)- X{z, t'-\) = U[,{Y{z, t')) , Vt' = t^^ + |, . . . , W - I 

Y{z,t + \)-Y{z,t-\)=\nq-Vl{X{z,t)) , Vt = t„,in, . . . , Wx - 1 

(4-5) 

2. There exists a point in C, which we call oo E C, such that X(z, tmin) has a simple 
pole at z oo, and we have: 

1 

F(^,t,,i,--)~^-— -. (4-6) 

3. There exists points Q G C, such that X((^j,tjnax) is an eigenvalue of Mf^^^: 

X iCu t^,^) = hi{iy), (4-7) 

and i^(-2, tmax — I) has simple poles at the points Q and behaves like: 

Y{z, Wx -h- ^, ^. , , . (4-8) 

^ X{z,trai,^) - hi{u) 
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4- Define for t = tmin, ■ ■ ■ , 



Y{z,t- ^)dX{z,t) 



W{x, t) = Res I^T^T^' (4-9) 

z^oo X-X{z,t) 

and call it "the resolvent" of the matrix Mf. We require that Vt, W{x,t) is 
analytical in a vicinity of x ^ oo, and behaves like: 

N 

W{x,t) - — , (4-10) 

x— >oo X 

and W{x, t) can be analytically continued to 

C\V{t) (4-11) 
where we recall that T>{t) C R zs a compact region ofM.. 

5. Typically, W{x,t) may have branchcuts or isolated singularities, like poles or log 
singularities. The set of points at which W{x,t) is not analytical is called the 
support: 

supp(t) = {x , H^(x, t) not analytical }. (4-12) 

The interior o/supp(t) is called "the liquid region" (it contains the cuts, it ex- 
cludes the isolated singularities) . We have supp(t) C T){t), and 

supp {t) C V{t) (4-13) 

IfT>{t) — Ui[ai{t),bi{t)], we require that^i: supp {t) fl [ai{t) , bi{t)] is connected. 

6. If the domain T>{t) at time t is a disconnected union V^t) = l-i^=i[ai{t),bi{t)], we 
require that \/ i — 1, . . . , m^; 

-/ W{x,t)dx^ni{t), (4-14) 



where the integration contour surrounds the interval [ai{t) ,bi{t)\ in the X{z,t) 
plane, in the clockwise direction. 

Finding functions satisfying all those requirements for a general domain, with gen- 
eral weights a{t'), (3{t') is a difficult problem. But for not too complicated domains and 
weights, some simplifications may occur, and we will see many examples of explicit so- 
lutions below. 

From now on, let us assume that we have found the functions X and Y satisfying 
all the requirements. Once we have found a solution to this problem, i.e. found the 
functions X{z,t) and Y{z,t'), we define the spectral curves: 
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Definition 4.1 The spectral curve at time t is the pair of functions: 

5, = (X(.,t),F(.,t-^)). (4-15) 

Remark 4.1 Because of eq. ()4-5p . the following 2-forms in T*C A T*C, restricted to the 
spectral curve, are equal: 

1 . , 1 

dX{z,t) AdY{z,t - -) = dX{z,t) AdY{z,t + -) 

= dX{z,t + l) AdY{z,t+-) 
(4-16) ^ 

and therefore, the spectral curves St and St+i are symplecticaly equivalent: 

5t = 5t+i. (4-17) 

4.3 Symplectic invariants and topological expansion 

The symplectic invariants Fg{S) were introduced in [2?]. To any spectral curve S, 
one can associate, by simple algebraic computations, an infinite sequence of complex 
numbers Fg{S), g = 0,1,2,3, . . .. We recall their definition in appendix [Bl 

One of their main properties, is that if two spectral curves S and S are symplecti- 
cally equivalent, then we have Fg{S) = Fg{S). 

In our case, because of eq. ( 14-17p . we have: 

Fg{St) = Fg{St+i), (4-18) 

and therefore Fg{St) is independent of t, the Fg(iSt)'s are conserved quantities. 

It was proved in [33J, for any chain of matrices, and here we apply it to our case, 
that: 

Theorem 4.1 

oo 

lnZ = 5^F,(5,) (4-19) 

9=0 

where the right hand side is independent oft. 

This theorem holds order by order in some appropriate formal large parameter 
expansion. We will see examples below, where the formal parameter can be the size of 
the system, or Ing, or a, ...etc. 
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4.3.1 Arctic circle 

For most interesting applications, there is some "large parameter Q" in our problem 
(typically the size of the domain T>, ot Q = 1/ Inq, or sometimes other parameters), 
such that the spectral curve scales like Q, typically: 

S = QS^ + o{Q) (4-20) 

(where we write XS = (x, Xy) for a spectral curve S = {x,y)). 

The large Q spectral curve S^o was already computed in many works and in par- 
ticular by Kenyon-Okounkov-Sheffield |l8j who found the limit shape of the liquid 
region. 

From the homogeneity property of F^'s (see [2^ and appendix [B]) we have Fg{S) = 
Q2-2g ir^(Q-ij5)^ i.e. we find a large Q expansion: 

oo 

lnZ = J2Q'''' FgiQ''^) (4-21) 

9=0 

Such an expansion is not very useful if Fg{Q^^ S) depends on Q. 

In fact, in many examples related to TASEP and plane partitions, we find that the 
spectral curve S depends on Q, but up to a symplectic transformation we have 
miraculously (this happens for instance in the matrix model considered in j26]): 

S = Soo modulo symplectomorphisms. (4-22) 

which implies: 

Fg{Q-'S) = Fg{S^), (4-23) 

and in that case, Fg{Q~^ S) = Fg{Soo) is independent of Q. 

This miracle is deeply related to the structure of the self-avoiding particles model 
partition function, and with the so called "arctic circle phenomenon", i.e. the fact 
that the system freezes beyond a certain size |12]. This can also be related to the fact 
that we have some arbitrariness in choosing the potentials Vt, and we could choose 
some potentials Vt which depend on Q in an appropriate way such that the spectral 
curve S would not depend on Q. Although it is doable in theory, finding the 
corresponding Vt's seems horrendous. A rather explicit example of this phenomenon 
was discussed in [26]. 

Only in the case where we have this "arctic circle phenomenon", we have: 

oo 

lnZ = Y,Q'~'' F'giS^), (4-24) 

9=0 
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where the spectral curve Soo is the curve derived from the Harnack curve of Kenyon- 
Okounkov-Sheffield ^8] . 

In that case, we can compute the large Q expansion of our self-avoiding particles 
model model, to all orders in Q, not only the large Q leading order limit. 

4.4 Reduced matrix model 

For a fully general domain V, there are as many equations eq. (14-50 to solve as the size 
of the domain, but, when the domain has very few defects, many of those equations 
simplify considerably, and we may consider a reduced problem. 

Notice, that everytime X>(t) = 0, we may choose Vt(x) = for all x, and thus 
V^' = 0. Therefore, it is possible to simplify dramatically the equations determining 
the spectral curve. 

Notice also, that there can be many times t at which 7^ 0, but since particles 
can only follow lines with slopes ±|, many places (the shadow of V) are never visited. 
In other words, we may replace the defects T> by the minimal defect Vq of T>, and since 
in general T>q is a very small subset of V, this allows to simplify the problem. 

Only the minimal defect T>Q{t) needs to be specified. 

Let the times at which Vo{t) 7^ be called: 

{^0) Ti, . . . , Tfc} C {tmiiu • • • ; ^max}; (^-25) 

among those T^'s we include the extremities tmin and tmax^ 

^0 = ^min ! ^fc = ^max- (4-26) 

The spectral curve equations eq. fl4-5p can be rewritten: 

X{z,T,+,)-X{z,T,) = E^=;;4 f/;,(f(^,T,+i - i) + (f -T,+i + i)lng) 

Y{z,Ti+, - i) - Y{z,T, - i) = (T,+i - T,) Ing - V^^{X{z,T,)) 

(4-27) 
We thus define: 

X,(z) = X{z, Ti) , Y,{z) = Y{z, Ti - ^) (4-28) 

and the potentials 

Ti-- 

U[{y)= J2 Ul,{y + {t'-T. + l)lnq) (4-29) 
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Figure 11: The left figure represents the domain T> and its defect T> in red. The 
right figure represents the minimal defect T>q in red. If there are k + 1 times Tq, . . . , 
containing red lozenges, then the reduced matrix model is a 2/c + 1 matrix model. 
The potentials of the reduced model must be such that e"^^* = at the center of red 
lozenges, and e"^^* = 1 at the blue ones. The values of e"^^* elsewhere are arbitrary. 



i.e. 

^-UM = Yl (a(t') e-f + (3{t') qU''~^^+\) e^) . (4-30) 

The loop equations eq. fl4-27p are equivalent to: 

J x,+i(^) - x,{z) = f/;+i(y;+i(z)) , = 0, . . . , - 1 

1 f,+i(^) - Y,{z) = (T,+i - T,) Ing - V^^iM^)) , = 0, . . . , A; - 1 

(4-31) 

One recognizes that these equations are exactly the equations of the spectral curve of 
another matrix model, which is a chain of matrices, but with only 2k + 1 matrices 
instead of 2(tmax — ^min) — 1, and we get the following result: 

Theorem 4.2 The matrix model partition function can be rewritten: 




k-1 



JJ^ Q-T^VTiiMi) ^(Ti+i-Ti) Tr Mi ^-TrUi(Ri) ^Tr Ri (Mi- Mi_i) 

i=0 1=1 1=1 

(4 - 32) 

i.e. a chain o/2/c + l matrices, instead of a chain of 2{tuiax — tmin) — 1 as we had before. 
Notice that knowing Xi, we can reconstruct X{z,t) for all times t e [Ti,Ti+i]: 

X{z,t)^M^)+ E U^,{Y{z,T,+,-^) + {t'-Ti+^ + ^)lnq), (4-33) 

and iffe [Ti + |,Ti+i-i] 

Y{z, t') = Y,+,{z) + (T,+i - ^ - t') In g. (4-34) 

5 Liquid region 

Our spectral curve is thus given by a collection of functions X{z,t) and Y{z,t + |), 
defined for all integer times t. We assume that the potentials Vt are non-zero only for 
times t e {Ti, T2, . . . , Tjt}, and the domain at time Tj is the union of intervals: 

V{Ti)^U[aij,k,j], (5-1) 

with filling fractions riij, i.e. we require that there arc riij particles in [(ii,j,bij]. 
This implies that the potentials V{x,Ti) are of the form: 

nii 

V'{x,T,) = f'{x,Ti) + J^^p{x- a„) - ij{x - 6,,,), (5-2) 

where f'{x, Ti) is some analytical function, and ip = T'/F is the digamma function. We 
used the well known property of ip that: 

n ^ 

: = ip{x) — ip{x — n). (5-3) 

i=l 

5.1 Interpolation to real times 

It is in fact better to enlarge these definitions to all times (not necessarily integers or 
half-integers) by interpolation: 

V t' e]r,_i, n , Y{z, t') = Y{z, - 1) + (T, - ^ - t') In q (5-4) 
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and we write: 

l/(z,t) = e^(^'*) (5-5) 



I.e.: 



V t' G]r,_i,r,], y{z,t') = yiz,T, - ^) g^-^*' = g^-^*' (5-6) 



The function y{z,t') is discontinuous at time Tf. 
Similarly we define for t G]Tj„i, Tj]: 



(5-7) 



If is constant over some interval, then we can find explicitly the interpolating 
function: 

t) ^ X(.,T.) + + ,,(-^ -1. _ (5.9) 

where 'ipq{x) is defined as (see appendix Rl) : 



If g = 1, we have: 



A-(..<)^A-(.,T.) + 421|^. (5-11) 



Notice that this expression depends linearly on t. 



5.2 Densities of particles 



Consider the density of particles at time t (in section 12.2.21 we have seen that these 
densities are related to the profile of the plane partitions): 



N 



p{x, t) = J2 (^(^ - W)) = ( Tr 6{x - Mt)) . (5-12) 



i=l 



We also consider their Stieltjes transforms, i.e. the resolvents: 
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Let us assume that our model has a formal large Q expansion in some parameter Q, it 
follows from the general solution of the chain of matrices [33], that one can write: 

oo oo 

W{x, t)=Y^ Q^-^' W^'\x, t) , p{x, t) = J2 Q^'^' P^'\^^ (5-14) 

5=0 9=0 

where the first term W^^\x,t) = W{x,t) is given by eq. f l4-9p : 

W^°\x,t) = W{x,t) = Res —-^ ¥ (5-15) 

z^oo X-X{z,t) 

and the density is the discontinuity of W{x,t) across supp(t): 

W{x-iO,t)-Wix + iO,t) 

^— = P{x, t). (5-16) 

All the higher corrections W^^\x,t) ioi g > 1 were also computed in [53], and they 
are the correlators Wi^^ of the symplectic invariants (see appendix [B] or [27]) of the 
spectral curve St- In other words, if we know the spectral curve St at time t, there is a 
simple recursive algebraic algorithm to compute all the corrections, of any correlation 
function, to all orders in Q. We shall not detail this here. 

Remark 5.1 Let us mention, that if the arctic circle property does not hold (this may be 
the case for complicated domains), then the spectral curve St may depend on Q, and thus all 
the W^^^ depend on Q. In other words, eq. (j5-14p is true as an equality of formal series, but 
it is not really a large Q expansion. In particular, W^^^ is not necessarily the large Q limit. 
When the arctic circle property holds, then, the W^^^ are independent of Q, and we really 
have a large Q expansion. We will see examples where this holds, below in section [71 

Let us now study the leading term. 

W{x,t) can have various sorts of singularities, it can have isolated singularities, 
like poles, or log singularities, and it can have branchcuts. In terms of the density 
p{x,t), the poles of W{x,t) correspond to 5-Dirac distributions, and cuts correspond 
to smooth positive densities. Log singularities of W{x,t) correspond to discontinuous 
jumps of the density. 

Notice that eq. fl4-10p implies that p{x, t) is normalized, with total weight A^: 

p{x,t)dx = N, (5-17) 

but if we exclude the isolated singularities, i.e. if we integrate only in the liquid region 
(the interior of the support), we have: 

p{x,t)dx<N. (5-18) 

supp(i) 
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Notice also that if T>{t) is disconnected, i.e. a union of intervals '^'(t) = 
U^\[aj(t), we have required that the filling fraction is fixed: 

/ p{x,t)dx = ni{t). (5-19) 
5.3 The Envelope of the Uquid region 

The liquid region, is the interior of the support of the density, i.e. it excludes all the 
isolated singularities of W{x,t), it contains only the cuts. 

The complement of the liquid region in V, is called the "solid" region, it contains 
the isolated singularities of W{x,t), and possibly their accumulation points. 

The boundary of the liquid region is called the envelope of the liquid region. It is 
given by the branchpoints. 




t 



Figure 12: Given a domain V, we compute its spectral curve St for all times t. From 
the spectral curves we compute the densities p{x,t), whose supports are contained in 
'D{t). The interior of the support is called the "liquid region". The liquid region is 
delimited by the envelope, which is a curve Xc{t) contained in T>. The points on the 
envelope are branchpoints of the functions X{z,t), i.e. zeroes of dX{z,t) We shall see 
below that the envelope has some special properties of tangency to the domain T>. 



The branchpoints Zc{t) at time t, are zeroes of the differential form dX, i.e. they 
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are the points where the density has a vertical tangent: 

dX{zc{t),t) = 



(5-20) 



where d is the differential with respect to z. Not all branchpoints are boundaries of 
the liquid region, let us consider only those which are on the envelope. The boundary 
of the liquid region is thus at Xc{t) such that: 

X,{t) = X{z,{t),t). (5-21) 

By definition, the envelope Xc{t) is a curve contained in V. Let us study some of its 
properties. 

5.3.1 Tangency 

Theorem 5.1 // a point of the envelope has a tangent of slope ±1/2, then this point 
must be on the boundary of the shadow ofV. 











/ \ ^'■^ 




Fi{ 


2;ure 13: A point of the envelope at which the tangent has a slope ±1/2, 


is necessarily 


on 


the boundary of the shadow of P, this means that there is a point 


of the defect 


V 


on the tangent. 





proof: 

Let t G]Tj_i, Tj] be such that the point (t, Xj(t)) is a point of the envelope at which 
the tangent has slope ±1/2 (resp. —1/2). Let a = Zc{t), i.e. we have Xc{t) = X{a,t). 
Since on the envelope we have dX = 0, we see that at all times t we have: 

dX,{t) _ dX{z,t) 
dt ~ di 



(5-22) 

z=z„(t) 
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and therefore we must have: 

dX{a,t) 1 



(resp. - -), (5-23) 



dt 2 ' ' 2 

and, from eq. (15-81) . this imphes that 

yi{a) = (resp.?/i(a) = oo). (5-24) 

Then, if yi{a) = (resp. oo), we see from eq. (15-81) . that Wt G]Tj_i,Tj]: 

X{a,t) = Xi{a) + (resp.Xi(a) - ^^). (5-25) 

and in particular at t = t, we have: 

X,(t) = X,{a) + (resp. X,(«) - ^^), (5-26) 

which means that the tangent to the envelope at point i, passes through the point 
a = X{a,Ti) at time Tf. 

a = X{a,Ti). (5-27) 

Then, notice that yi{a) = e-*^'-"'-^*"^) can be equal to or oo, only if Y{z, Ti — |) has 
a singularity at z = a. Because of eq. (14-3 ip . we see that, singularities of Y{z,Ti — |) 
can occur only if V'{a, Ti) has a singularity. 

Remember that V'{x,t) has no singularity on T'(t), it has singularities only at the 
defects Vit), and thus, this means that the point (a, Tj) belongs to the defect V. But 
we know that the point {Xc(t),t) belongs to V, and the slope between the points (a, Tj) 
and (Xc(f),t) is equal to ±i, i.e. the point {Xc(t),t) is in the shadow of the defect 
{a,Ti}. 

□ 

5.3.2 Local Convexity 

In the classical case, where q*a{t) / [3{t) is constant, we have a convexity property: 
Theorem 5.2 Ifq*a/(3 is constant in time, the liquid region is locally convex. 
proof: 

The density at time t must have a real support, starting at a branchpoint, i.e. at 
Xc{t). Consider that, up to a reparametrization of z, and up to trivial translations in 
t and X, we can write localy, near the branchpoint point Xc(0) we write: 

X{z, t)=X, + Xtt + ^X,,,z2 + x.^tzt + ... (5-28) 
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and up to a rescaling of we assume Xz^z = 1. Notice from eq. (15-111) . that the 
assumption q*a/(3 constant guarantees that X{z,t) is hnear in t, there is no term. 
The branchpoint at time t is at Zc{t) given by dX{z,t)/dz = 0, i.e.: 

z,{t) = -X,^t t + ... (5-29) 

i.e., the envelope Xc(t) = X{zcit),t) is, up to order 2: 

X,{t)=X, + Xtt~^e + Oit^) (5-30) 

Our assumption on the reahty of the supports, imphes that Xc{t) G M for all t, and 
therefore Xt and X^^ must be real quantities. 
The point z such that X{z,t) = X{z,t) is: 

z = 2z^{t) -z + ... (5-31) 

The function y{z) is then given by solving eq. (15-81) : 

. . _ J-T,.^l <t') 1-2{X, + Xz,z + ...) 
' ' ~ ^ m l + 2iX, + Xz,tz + ...) 

- ^ m 1 + 2X, 1 - 4X2 +■■■) (5-32) 

and therefore 

y{z,t') ^ 8X,,t 



^2 {z - z,{t)) + ... (5-33) 
-t 

and the density p(x, t) = 2^ In with x = X{z, t) is, up to order 2: 



y{-z,t') 1-4X1 



This implies that X^^ i^dt) — a;) > in the liquid region. If the liquid region is above 
Xc(t), we must have Xf ^ < 0, and if the liquid region is below Xc(t), we must have 
X^^ > 0. In all cases, the liquid region is locally convex. □ 

Remark 5.2 If we are not in the classical case, i.e. if q^a/fi is not constant, then X{z,t) 
is not linear in t, there maybe a second derivative Xt^t 7^ 0, and the envelope is locally: 

X,{t) =X, + Xtt + ~ f + Oif) (5-35) 

and although X^^ has a constant sign, Xt^t — X^^ can have any sign, and the envelope is not 
necessarily convex. 



38 



Figure 14: The envelope has a cusp at t = Tj. 



5.3.3 Cusps 

Notice that Y{z, t) is discontinuous at t = Tj, i.e. there is a left value and a right value. 
This implies that dX^/dt has a right tangent and left tangent, and thus generically, 
the envelope has cusps at t = Tj. 

5.3.4 Genus and holes 

Condition 6 of section 14. 2[ i.e. eq. fl4-14p tells us that the spectral curve must have at 
least as many non trivial cycles as the number of holes in our domain, i.e. the genus 
of our spectral curve is at least the number of holes, and generically, it coincides with 
the number of holes. 

In particular, if we consider a simply connected domain T) with no holes, we can 
expect to have genus 0, which means that the functions X{z,t) and Y{z,t) are ana- 
lytical functions of a complex variable z. The Riemann surface in which z lives can be 
chosen domain of the complex plane. 

6 Asymptotic regimes 

Of course, our model depends on so many parameters (shape of the domain, coefficients 
a{t'), Pit'), potentials Vj, ...) that it is almost impossible to classify all possible asymp- 
totic regimes. However, a few asymptotic regimes are more relevant for applications in 
statistical physics or algebraic geometry. 
We classify them into 2 kinds: 

• macroscopic asymptotic regimes, which describe the behavior of the partition 
function, i.e. the statistics of our self-avoiding particles model at the scale of the size 
of the domain. Typically, we shall consider asymptotic behaviors for large size, of for 
small Ing. 

• microscopic asymptotic regimes, which describe the statistics of our self-avoiding 
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particles model particles in a very small region of the domain, typically, the behavior 
in the bulk, or near the edges, especially near special points, like near the envelope, or 
near cusps of the envelope. 

Let us comment a few of them, and let us emphasize that our method, works for 
all possible asymptotic regimes, and it gives not only the leading order asymptotics, it 
gives the full asymptotic expansion to all orders. 

6.1 Classical case q=l, and large size asymptotics 

Classical means that we choose q = 1- Here, we shall assume that the weights a{t') = 
P{t') = 1 are constant, although the general case is doable by the same methods. 

Consider a domain V whose size is T = tmax — ^min- Assume that there are k defects, 
and k does not depend on T: 



The defects are at times Tj = Tx;, where Xj are independent of T, and the defect at 
time Ti is the union of rrii intervals: 



k = 0(1). 



(6-1) 



(6-2) 



where again aij and bi j don't depend on T. 

Assume also that the number of particles scales like T: 



N = J\fT, 



(6-3) 



and the number of particles rij .,- in each interval scales like T: 



(6-4) 



and of course Vz 




(6-5) 
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See fig [13 



6.1.1 The spectral curve 



The spectral curve equations eq. (14-3 ip are: 




(6-6) 



where we choose the potentials Vi according to eq. (13-141) : 



Vl{X) = f\X) + 5^^(X - Ta,,,) - V^(X - Th,) 



(6-7) 
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Figure 15: A polygonal domain with k defect at times Tj = Trj. The envelope of the 
spectral curve, is the algebraic curve of smallest degree, tangent to all boundaries. 



where /' has no singularity, and ip = T'/T is the digamma function. We used the well 
known property of ip that: 



n ^ 

: = ijjix) — iljix — n). (6-8) 



X — 2 
i=l 



Moreover, the last matrix of the chain = M(tniax) = diag(Aj) is not integrated 
upon. Its spectrum is ^'(tmax), i-e. it is delimited by the intervals [Takj,Tbkj]- Then, 
equation eq. fl4-8p reads: 



y ~ , M EE ^r—^ =J2^P{X,- Ta,,) - ^PiX, - Th,) (6-9) 

Xfc^e.v.ofMfe ^ . Xk — I ^ 



6.1.2 Rescaling 

Since our model is defined as a formal power series, and we are looking for the large T 
expansion, we should try to find the spectral curve in the large T expansion. In that 
purpose we rescale the variables: 

Xi = Txi , Yi = \ny,, (6-10) 
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and we use Stirling's asymptotic formula for the F and ip functions, see appendix |X1 

1 °° B 

^{x)r^lnx---J2^, (6-11) 

n=l 



where Bn are the Bernouilli numbers. 
That gives: 

rrii , 

^ = J] ^IJL^ (1 + regular + 0(1/T)) (6-12) 

where the 0(1/T) term is, to each power of T~^, a rational function of Xi with poles at 
ttij or 6, J, plus possibly an arbitrary analytical function with no singularity. We also 
have: 

nik 

Vk-^U ~ y (1 + regular + 0(1/T)) (6-13) 
j=i 

where the 0(1/T) term is, to each power of T~^, a rational function of with poles 
at akj or bkj, plus possibly an arbitrary analytical function with no singularity. 
We also have from eq. (16-61) : 

Xi - Xi_i = — — 6-14 

2 l + Vi 

i.e. 



_ n - r,_i - 2{xi - Xi_i) 

Ti - Ti^i + 2{Xi - Xi-i) 

6.1.3 Singularities 



(6-15) 



We see from eq. (16-121) .and eq. (16-131) . that if Xi = aij, the ratio yi/uij^i has a zero, and 
thus either Ui has a zero, or yi^i has a pole, or both. 

In fact, since Ui and Ui+i are multivalued function of Xj, there might be several 
points on the spectral curve such that Xj = aij. In that case, it may happen that 
Hi = corresponds to Xi = in one sheet, and = oo corresponds to Xi = Qij 
in another sheet. Or if both yi has a zero, and has a pole at the same point, this 
means that Xi — aij has a double zero, i.e. dxi has a zero, and this means that we are 
at a branchpoint, i.e. the point where two sheets meet. 

Therefore, we shall assume that generically both possibilities occur, and we conclude 
that: 

• Hi has a zero when Xi = ajj or = and has a pole when Xi = bij or 

Xi—i Cli—ij. 

Since the Xj's and ?/i's have only meromorphic singularities, it is natural to look for 
an algebraic solution, i.e. such that Xi and i/i are meromorphic functions on an algebraic 
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curve L. Indeed, our freedom to choose the potential Vt can be used to ehminate all 
the singularities except those described above. 

The problem then consists in finding algebraic functions Xj, such that: 

-^i -^i-l — 2 

• Di has a simple zero when Xi = ai,j or = and has a simple pole when 

• the density measures ^ In ^^^j dxi{z) (where = Xi{z) are the two points 
on each side of the support) must have their supports in V (in particular the supports 
are real). Those reality conditions for an algebraic curve, are closely related to the 
conditions which defines the Harnack curve in Kenyon-Okounkov-Sheffield [48] . 

6.1.4 The large size limit: Harnack curve 

So, we have to find a spectral curve, and its envelope, satisfying many equations (which 
give the poles and zeroes), as well as many reality conditions. Since we have only mero- 
morphic types of singularities (poles), it is natural to look for real algebraic spectral 
curves. 

It is also natural to look for minimal degree algebraic curves, i.e. having just the 
number of poles and zeroes implied by our equations and no other poles. In fact 
introducing other poles would most probably break the condition on the homology of 
steepest descent paths. Also we can again use our freedom to choose the potential Vt 
in order to reduce the degree, so that we eliminate all the singularities except those 
described above. 

An envelope and spectral curve having all the required properties can be found 
from the work of Kenyon-Okounkov-Sheffield |18], and they proved that their envelope 
is indeed the limit shape of the liquid domain to leading order at large T. 

Their spectral curve is a Harnack curve, and the envelope Xc(r) is an algebraic 
curve given by a polynomial equation 

P(xe,r) = (6-16) 

where P is a real polynomial, which has some very special properties. 

In particular, it satisfies all the properties we want for our envelope, and also, being 
a maximal Harnack curve means that the area enclosed by the envelope is maximal. 
Another property, is that the genus of P, is exactely the number of holes of the domain 
in the envelope. 
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6.1.5 Recipe for an algebraic spectral curve 

Therefore, let us make the assumption that the spectral curve is the algebraic curve of 
smallest degree satisfying all the constraints. We shall discuss the consistency of that 
assumption in section [5.1.71 below. 

The problem we have to solve in order to find the spectral curve is then: 

• Find a Riemann surface whose genus is equal to the number of holes in V. 

• Find 2 meromorphic functions u{z) and v{z), with only simple poles, and let us 
denote the set of poles as cxDjj, i = 0, . . . , /c, j = 1, . . . , rrii. 



We denote 
and Vi: 

Xi{z) = U{z) + TiV{z). 

and we require that Vi,j (cusp condition): 

Res Xi{z) = 

The meromorphic functions u{z) and v{z) must be such that Vi, j: 

1 



if z < A; 3z such that Xi(z) = aij and viz) = — 



if i > 3 2; such that Xi{z) = aij and v{z) = +-, 

ifi < k 3z such that Xi{z) = bij and v{z) = +-, 

ifi > 3z such that Xi{z) = bij and v{z) = —-. 
and we have the filling fractions: 

7^ [ ^dy{z)=Ar,, 

where 

, , l-2v(z) 



l + 2v(z)' 



6-17) 



6-18) 



6-19) 



6-20) 

6-21) 
6-22) 
6-23) 

6-24) 
6-25) 



and where we assume that the number of particles in [ciij, bij] is riij = T Mij. 
Explicit examples for given domains (hexagon, cardioid) are treated in section [71 
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6.1.6 Envelope of the liquid region 

The envelope of the liquid region corresponding to the spectral curve above, is obtained 
as follows: 

the branchpoints Zc{t) are solutions of dx{zc, r) = = du{zc) +Tdv{zc), and in fact, 
it is much easier to compute function of Zc, namely: 

Then, one computes Xc(t) = x{Zc{t),t), and again, it is easier to parametrize this 
equation by Zc, that gives: 

f ,M^.) ■ (6-27) 

We see that the envelope is an algebraic curve, and from section 15.3. ll we know that it 
is an algebraic curve tangent to the boundary of the shadow of T>, with cusps at r = r^, 
and whose genus is the same as the number of holes of T>. 

Again, explicit examples of envelopes are given in section [71 For example, the 
envelope of an hexagonal domain is the ellipse tangent to all sides of the hexagon, see 

fig [13 

Recovering the spectral curve from the envelope: 

Assume that we know the envelope Xc{tc), i.e. the algebraic curve tangent to the 
boundaries of T>. Let us explain how to recover the full spectral curve x{z, r). 

Given the equation of the envelope Xc{tc) (which is a multivalued algebraic func- 
tion), one can choose locally Z Tq clS db local parameter. One finds that the spectral 

curve is (at least in the domain of r where Tc can be chosen as a local parameter): 

x{z,t) = Xc{z) + {t — z) v{z), 

I \ \-2v(z) 1 / \ dxAz) (6-28) 

Therefore, knowing the envelope allows to recover the full spectral curve, i.e. the 
functions x{z^t) and y{z). 

6.1.7 Large size Asymptotic expansion 

Now, suppose that we have found those functions, and that we have indeed found 
the correct spectral curve. The spectral curve is then the pair of functions St = 
{X{z,t),Y{z,t)), i.e., up to a symplectic transformation: 

S=iTuiz),\niy{z))), (6-29) 
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and we notice that u{z) and y{z) do not depend on T, therefore we have Fg{Tu, Iny) = 
T^-^>^Fg{uMy). 

And thus, from theorem I4.H we have the full large T asymptotic expansion: 
Conjecture 6.1 InZ has the large size T expansion 

oo 

InZ ~ ^T2-2sF^(^i,ln2/). (6-30) 

5=0 

where the meromorphic functions u{z) and y{z) are computed by solving the require- 
ments of section \6.1.5\ and Fg is the g^^ symplectic invariant defined in l2^ . 

sketch of a possible proof: 

A hint to that conjecture, is that the leading large T densities are governed by 
Fo(m, Iny), and it can be seen that this coincides with the limit shape found by Kenyon- 
Okounkov- Sheffield [i8] . 

In order to prove this conjecture, relying on the work of [33], we only have to prove 
that we have indeed found the correct spectral curve, and that our guess (that u and 
V are algebraic functions of smallest possible degree) is correct. 

In principle, this means proving that the integration paths in our self-avoiding 
particles matrix model, are indeed the steepest descent paths for our potentials, but 
this seems too difficult. For matrix models with polynomial potentials, this is usually 
proved by the Riemann-Hilbert method of Deift & co [17] . 

Another possible proof, is to prove this order by order in some formal parameter, 
especially if the model tends to a Gaussian matrix integral in the small parameter limit. 

One suggestion is to notice that in our matrix model, the potentials Vt may depend 
on T, and our spectral curve can be expected to depend on T, and somehow, the 
Harnack curve gives only the leading term: 

S{T) = 5oo + 0(1/T) (6-31) 

where Soo is given by the Harnack curve of |48j . Indeed, we have not really taken into 
account the 0(1/T) term in eq. fl6-12p and eq. (16-131) . 

However, we have some freedom in the choice of VJ, and we can choose any Vt 
provided that the spectral curve satisfies eq. (16-121) and eq. (I6-I3p . and in particular, 
we may choose the spectral curve S^o, which does satisfy eq. (16-121) and eq. (16-131) . 
with the 0(1/T) term vanishing. In other words, we use the freedom of choosing the 
potentials Vt, and we define Vt from the spectral curve through eq. (I6-I2p . instead of 
the contrary. 
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Somehow, we go backwards, we first construct the spectral curve iSoo, and then we 
construct the potentials Vt which correspond to it. 

This very special choice of Vt guarantees that Soo is the spectral curve of our model, 
and then it is independent of T □. 

6.2 Quantum case 

Now we consider q ^ 1. We also assume for simplicity that the weights a{t') = (3{t') = 1 
are constant. 

The potential UiiY) appearing in eq. fl4-30p is: 

(6 - 32) 

where g{y) = YYjLii^ " Q'' /v) is the g— product (it is a quantum deformation of the 
F-fimction, see appendix |A]). And thus: 

Um = ^1^^ +^,(-e-^) -^,(-e-^g(^-^-^)) (6-33) 

where ipq{x) = xg'{x)/g{x) = J2'jLiQ''/{^ ~ Q'"')- 
Thus, eq. (14-311) gives: 



v V _ Ti^l 

■^i — ^i-i — 

and at intermediate times Tj_i <t<Ti: 



E y (6-34) 



Ti-t-l 



We choose the same domain T) as in the classical case, i.e. a domain which scales 

with a factor T = tmax — ^min- 

Several possibilities may occur: 

• The regime T ^ j^, is more or less the same as g = 1, which we have studied in 
the previous section, i.e. there is a liquid region of typical size T. 

• In the regime T 3> there is a liquid region of typical size 1/| lng|, and most 
of the domain D is in a frozen phase. 

• In the intermediate regime ~ 0(1), the liquid phase if of typical size T ~ 
This is the most interesting regime. 
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6.3 Case - 0(1) 

We consider the regime where T is large and Ing is small, and g-^ ~ 0(1)- We define: 

= q (6-36) 
Then we shall repeat most of the steps of the classical case q = 1- First we rescale: 
x, = q^^ , yi = e^^ , T = t/T. (6-37) 

6.3.1 Equation of the spectral curve 

The equations eq. (14-3 ip read: 

"^i bi ■ 

yi±l ^ ^in+i-n) TT " ^ ^ ^^^^^^^ ^ O(lng)) (6-38) 



n ^' ~ (1 + regular + O(lng)) (6-39) 



Vk 



,=1 - 



And to leading order at small Ing we have: 

^ = q-"^ ^i±qr:il (1 + o(lng))) (6-40) 

Notice that those equations imply only meromorphic singularities for yi and Xj, 
and again, it is natural to look for an algebraic spectral curve. Notice that if yi has 
a zero (resp. a pole) at z, then we have Xi+i{z) = q^^'^^+'^^'^^^Xi{z) (resp. Xi_i{z) = 
qi■^^-l~■^^)x^{z)). 

6.3.2 Recipe for an algebraic spectral curve 

Therefore, let us make the assumption that our functions Xi{z) are algebraic of smallest 
possible degree satisfying all the constraints. The consistency of that assumption can 
be discussed like in section 16.1.71 for g = 1, and we discuss it again in section 16.3.41 
below. 

The problem we have to solve in order to find the spectral curve is then: 

• Find a Riemann surface whose genus is equal to the number of holes in V. 

• Find 2 meromorphic functions u{z) and v{z), with only simple poles, and let us 
denote the set of poles as ooij, i = 0, . . . , k, j = 1, . . . , rrii. 
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We denote 

x{z, r) = q^^\{z) + q-^/^v{z), (6-41) 

and Vi: 

Xi{z) = x{z, Ti) = q^^/'^u{z) + q-^'/^v{z). (6-42) 
and we require that Vi,j (cusp condition): 

Res Xi{z) = (6-43) 

The meromorphic functions u{z) and v{z) must be such that Vi,j (tangency 
conditions) : 



\ii<k 3z such that Xi{z) = q"^-^ and Xi+i{z) = q°-'-^ 2 ^ 

(6-44) 

if i > 3z suchthat Xi{z) = q°-^'^ and Xi_i{z) = q"-^-^^ 2 , 

(6-45) 

if i < /c 3 2; suchthat a;i(2;) = q and a;i+i(2:) = q 2 , 

(6-46) 

2 

(6-47) 



ifi>0 3z suchthat Xi{z) = q ''' and a;i_i(2;) = q ~ 



1 

2i7r 



and we have the filhng fractions: 

Yi{z)dX,{z)=nij (6-48) 
where Xi = q^' and Yi is obtained by solving eq. (16-341) . 

Exphcit examples for given domains (hexagon) are treated in section [71 
6.3.3 Envelope 

The envelope is given by the branchpoints Zc{t) solutions of dx{zc, t) = 0, and again, 
it is easier to find r as a function of Zc than the contrary. We have: 



The envelope Xc{t) is also better given in a parametric form with the parameter Zc as: 

(6-50) 



t _ _ dv(zc) 
y du(zc) 
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Since u and v are meromorphic, this implies that g* and q^" are related by a polynomial 
equation: 

P{q\q^-)=0. (6-51) 

Again, we claim that this polynomial is the same as the Harnack curve of Kenyon- 
Okounkov-Sheffield |18]. 

In other words, the plane curve Xc{q*) = q^''^^\ is an algebraic plane curve inscribed 
in the image of the domain V under the map (x, t) {q^, g*). 




X 

J ^ 



Figure 16: The envelope of the liquid region, is the smallest degree algebraic plane 
curve, tangent to the image of the hexagon under (x, t) (g^, g*). 



Recovering the spectral curve from the envelope: 

Assume that we know the envelope Xc{t), or, writing Xc = q^" and z = g*/^, assume 
that we know the plane curve Xc{z), i.e. an algebraic curve tangent to the boundaries 
of the image of V under {x,t) (g^',g*/^). Let us explain how to recover the full 
spectral curve x{z,t). 

Given the equation of the envelope Xc(g*^^) (which is a multivalued algebraic func- 
tion), one can choose locally z = g*/^ as a local parameter. 

One finds that the spectral curve is (at least in the domain of t where g* can be 
chosen local parameter): 

x{z, t) = ^ + x'^iz)^ + ^ (zx^iz) - z'x'Xz)) . (6-52) 

Therefore, knowing the envelope allows to recover the full spectral curve, i.e. the 
functions x(z, r) = g^(^'*) and also Yi[z) from eq. (16-341) . 
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6.3.4 Small Ing Asymptotic expansion 

Now, suppose that we have found those functions, and that we have indeed found 
the correct spectral curve. The spectral curve is then the pair of functions St = 
{X{z,t),Y{z,t)), i.e., up to a symplectic transformation: 

5 = (^j^ In {x{z, r)). In {y{z, r))^ , (6-53) 

and we notice that x{z,t) and y{z,T) do not depend on q, they depend only on q, 
therefore we have Fg(j^ Inx, Iny) = (lng)^^~^Fg(lna;, Iny). 

And thus, from theorem 14. H we have the full large T asymptotic expansion: 

Conjecture 6.2 InZ has the small Inq expansion 

oo 

\nZ ~ '^{\nqf3-'^Fg{lnx,lny). (6-54) 

9=0 

where the meromorphic functions x{z) and y{z) are computed by solving the require- 
ments of section \6. 3.2\ and Fg is the g^^ symplectic invariant defined in l2l^ . 

A possible proof could follow the same ideas which we discussed in section 16.1.71 for 
q=l. 

6.4 Microscopic asymptotics 

It is a truism to say that since our model is a matrix model, it has all the local universal 
behaviors of matrix models. 

6.4.1 Zoom near a point 

Let us choose a point (xq, to) anywhere in the plane (it can be in or outside the domain, 
or on the border). 

Let us rescale x and t with some scaling parameter s, with some exponents a and 
5, i.e. we write: 

x = xo + s° ^ , t = tQ + s^T. (6-55) 

Also, on the spectral curve, we choose a point zq such that X(zo,to) = a;o, and we 
rescale it by choosing a rescaled local parameter C,: 

z = Zq + s'~^ C,. (6-56) 

We rewrite the asymptotics of the functions X{z,t) and Y{z,t) in terms of those 
rescaled variables: 

X{z, t) = xo + s"e(C, r) + o(s°) , y (z, t)=yo + s^fj{(, r) + o(/). (6-57) 
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If the exponents a,/3,7,5 are well chosen, then the spectral curve iSzoom = 
(^(C) ^(C) ''"))) is a regular spectral curve (it has only branchpoints of square-root 
types). We call this curve, the blow up of the spectral curve {X,Y) near the point 
{xQ,tQ). In practice, finding the blow up is a rather trivial task, it merely consists in 
finding the first non- vanishing terms in the Taylor expansion near a point. 

From p7j, we have that the correlation functions Wn\xi, . . . , x„; t): 

^ Wjf^ (xi, . . . , t) = / tr tr . . . tr (6-58) 

^ \ Xi - Mt X2 - Mt Xn- Mt/ ^ 

behave at small s like: 

Wi^\x„ . . . , t) ~ ("+'3)-"" 4^)(ei, . . . , a; r) (1 + o(.°)), (6-59) 

where Un^ are the "symplectic invariants" correlation functions of [27] associated to 
the spectral curve iSzoom- 

In other words, since the correlation functions are symplectic invariant correlators 
of a spectral curve S, their local behavior is, to leading order, given by the symplectic 
invariant correlators of the blown up spectral curve. This theorem found in [27j is very 
easy to prove by recursion on n and g (see appedix iBj) . 

Therefore, it suffices to find the blown up spectral curve to characterize the leading 
behavior of the correlation functions near a point. 

6.4.2 Airy kernel near regular boundaries of the liquid region 

Near a regular point of the envelope Xdio), we have dX/dz = at t = and thus, 
we have a Taylor expansion in z of the form: 

X{z, t) ~ Xo(t) + {t- to)Xi{t) z + X2{t)z^ + . . . , (6-60) 

where each Xi{t) has a regular Taylor expansion in t — to- Let us rescale: 

z = sC. , t = to + ST. (6-61) 

We have to the first few orders in s: 

X{z,t)=Xo{to)+s (XoT^+s' (XirC + X2C')+0(s3). (6-62) 

We also have 

X{z, t)=Xo + s (XiCr) + 0{s^), (6-63) 

and therefore 

Y{z,t) r^Yo + sXiC + 0{s^). (6-64) 
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The blown up curve is thus: 



I yiC,r) = C 

which satisfies: 



(6-65) 



{y + T/2f = X + rV4. (6-66) 

The spectral curves such that Polg(y) = Polp(a;) where Polg is a polynomial of degree 
q and Polp is a polynomial of degree p, appear in the so-called (p, g)-minimal models, 
i.e. in the classification of conformal field theories [l9l [20]. It has central charge 
c = l — Q{p — qY/pq. Here, this is the (1, 2) minimal model, with central charge c = —2, 
which is well known to be generated by the Airy differential system Ai" = xAi, and 
the correlation functions are determinants of the Airy kernel. It is also well known to 
be related to the Tracy- Widom law of extreme eigenvalues statistics [70] . 

6.4.3 Pearcey kernel near cusps 

We have a cusp each time a pole disappears, generically a simple pole. Thus we have a 
Laurent expansion in z starting at z~^, and such that the residue of the pole vanishes 
at t = to, i-e.: 

X(z, t) ~ ^^"^"^""^^^ + xo(t) + u{t)z + ..., (6-67) 

and all the coefficients v{t), xo(t), u{t), . . . have a regular Taylor expansion near t = to- 
Let us rescale z = s(/u(to) and t = to + s'^T/u(tQ)v(to), to order s, and up to constant 
terms, we have: 



X{z,t)=xo + s {C + ^] + 0{s'), (6-68) 



and generically, Y behaves like: 



Y{z,t)=yo -- + .... (6-69) 

The blown up curve is thus: 

_ r x(C,r) = C+f 

<->Pearcey — S \ > 5 {K)-lU 

and the local behaviors of correlation functions, are the correlation functions of that 
universal curve. 

This is the spectral curve whose correlators are generated by the Pearcey kernel. 
We thus recover the well known Pearcey kernel behaviour fTT]. 
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6.4.4 Critical points 

If we consider a point on the boundary, such that more derivatives vanish, typically we 
find that the spectral curve behaves locally like: 

P 

y{x) ~ y{a) + (x — x{a)) « + . . . (6-71) 

It was shown in p7], that the correlation functions tend towards those of the (p, g) 
reduction of the K-P hierarchy, i.e. the (p, q) minimal model of central charge c = 
1 — Q{p — /pq. 

The model (3,2) of central charge c = is called "pure gravity", the model (5,2) 
of central charge c = —22/5 is called Lee- Yang, The model (4,3) of central charge 
c = ^ is called the Ising model,... Their correlators are generated by determinantal 
formulae from a kernel involving the ?/^-system associated to a non-linear equation of 
Gelfand-Dikii type (Painleve I is the Gelfand-Dikii equation for pure gravity (3,2)). 
Some details can be found in [201 E]- 

6.4.5 Other local behaviors 

Also, we expect that locally in the bulk of the liquid region, the behavior is given by 
the sine-kernel |TOl [6l |13], and near vertical boundaries of the liquid region, we have 
y{x) ~ y{a) + {x — a)^, i.e. we expect the (2, 1) model, described by the Hermit kernel 
of the "birth of a cut" (see [131 El])- And we also expect to find the "Bead model" 
limit in the tentacles of the amoeba when |Tlng| is large, see [12] . 

However, these cases are such that the blown up curve is not regular, and we cannot 
directly apply the method of [27|. 

6.4.6 Arbitrary local behaviors 

Then, one could easily invent some domains, for which a local blown up curve could 
be any spectral curve specified in advance, and by choosing sufficiently complicated 
domains one can obtain any limit law. 

The classifications of all possible laws is more or less the classification of singularities 
of spectral curves, and it is more or less the classification of spectral curves themselves. 

7 Examples 

In this section, we illustrate our method by applying it to several classical examples, 
which were already studied in the literature with other methods. Here, however, we 
have a method to obtain not only the large size limit, but also all corrections to all 
orders. 
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7.1 The hexagon 

Our domain V is the hexagon of figure WJ\ with slopes ±|. 




We choose tmax = —tmin = T. We choose the weighs a(t) = P(t) = 1 at all times. 
We write N = bT. 

There is no defect, therefore we have k = 1, the reduced matrix integral is a 2- matrix 
model with an external field: 



-2^hexagon — / dMf) 



dRi Q-T^^VoiMo) g-Tr!7i(i?i) gTrRi(Mi-Mo) 



(7-1) 



where Mi is the fixed following matrix: 

Mi = dmg{Ta + l,Ta + 2,...,Ta + N) = T diag(a + ^, a + 



For Vq, we can choose 



-Vo(X) 



1 



T{X)T{N + 1- X) 
7.1.1 The classical hexagon q = 1 



a + b). (7-2) 



(7-3) 



We apply the recipe of section I6.1.5[ 

First, the domain has no hole, and we look for an algebraic curve of genus 0, 
therefore it can be parametrized by a uniformizing variable in the complex plane z G C. 
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Let us look for 2 rational fractions u{z) and v{z), with 2 poles, and we write: 

x{z,t) = u{z) + Tv{z). (7-4) 

Up to a reparametrization of z, we choose the poles to be at z = and z = oo. 
Moreover, we require that at r = — 1 the pole at ^ = disappears, and that at r = +1 
the pole at 2; = oo disappears, this implies that u and v and x are of the form: 

x{z, r) = c + rr + 7 ^(1 - t)z + 5 (7-5) 

u{z) = c + 'J ^z + ^ , v{z) = r + J ^-z + ^ . (7-6) 

We define: 

Xo{z) = x{z, —1) = c — r + 2'yz , Xi{z) = x{z,l) = c + r + 2'y / z. (7-7) 
Then we find the coefficients r, c, 7 by solving the following system: 



3z 


such that 


Xoi 


> 


3z 


such that 


Xoi 


> 


3z 


such that 


Xii 


> 


3z 


such that 


Xii 





and f (z) = — 



6 and f (z) = i, 
Xi(z) = a and v{z) = ^, 
xi{z) = a + b and v{z) = 



(7- 



We have 3 unknowns and 4 equations, but one can check that there are only 3 inde- 
pendent equations, and the system has a solution. 
We easily find: 

a + b a(l + b) 9 , , n w, n / n 

c=^ , r = ^^ , 167'= (l-a)(l + «)&(& + 2). (7-9) 

Envelope 

The branchpoints Zc{t) are found from x'{zc, r) = 0, and we find that there are two 
branchpoints: 

^c(r) = (7-10) 

That gives: 

x^{t) = c + rr ±27^1 -r2, (7-11) 

or explicitely 



^ (^a + 6 + ra(l + 6) ± 7(1 - r^) (1 - a^) 6(6 + 2)) . (7-12) 



One can easily check that this is the equation of the ellipse tangent to all sides of the 
hexagon. See the ellipse in figJTTl 
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7.1.2 The quantum hexagon 

Now consider = 0(1), and write: 

q = / (7-13) 

and we shall define the q-numbers: 

W = (7-14) 

q 2 — qz 

Let us now apply the recipe of section I6.3.2[ 
We write: 

x{z, r) = q^/^u{z) + q'^/%{z) (7-15) 

where u and v are rational fractions with two poles. Up to a reparametrization of z, 
we choose the poles to be at 2; = and z = 00. Moreover, we require that at r = — 1 
the pole at 2; = disappears, and that at r = +1 the pole at z = 00 disappears, this 
implies that u and v and x are of the form: 

r) = {c + z) + Il±^ (r + ^), (7-16) 



We define: 



Xo{z) = x{z,—l) = c + z , Xi{z) = x{z,l) = r + -. (7-17) 



Then we find the coefficients r, c, d by solving the following system: 



3z 


such that 


Xo{z) 


= 1 


and 


Xi{z) 


= q-\ 


3z 


such that 


Xo{z) 


= q^ 


and 


Xi{z] 


1 = q^+\ 


3z 


such that 


xi{z) 


= q^ 


and 


Xo{z] 


) = q"-\ 


3z 


such that 


Xi{z) 


= q'^+^ 


and 


Xo( 


'z) = q-^+^+i 



(7-18) 



We have 3 unknowns and 4 equations, but one can check that there are only 3 inde- 
pendent equations, and the system has a solution. We find: 

q-i ^ qa+fe _ _ qb+i 1 + q«+b+i _ q^-i _ q^ 



q-i-q ' q^^-q 



(7-19) 



(q-i -q)2 



Envelope 



The branchpoints are at Zc{t) = ± a/c? [1 + r] / [1 — r] , and thus the envelope is: 
See some examples plotted in figure 
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Figure 18: The envelope of the hexagonal domain. The plots of eq. ( l7-2ip are for the 
hexagon a = 0.3, b = 2, and for values of q respectively q = 0.001, 0.1, 0.3, 0.8, 10, 1000. 
Notice that the liquid region is convex only for q close to 1. 



7.2 The cardiod 

Consider the classical case q = 1- The domain is the one represented in figure [191 We 
assume < a < 1, and 6 > 0. 

We find the spectral curve by applying the recipe of section 16.1.51 



where u and v are rational fractions with 3 poles. 

Up to a reparametrization of z, we choose the poles to be at 2; = —1, 0, 1. Moreover, 
we require that at r = — 1 the pole at 2; = — 1 disappears, at r = the pole at z = 
disappears, and at r = +1 the pole at 2; = +1 disappears. Moreover, since the domain 
has a symmetry r ^ — r, we choose u and v and x of the form: 



Consider 



x{z, t) = u{z) + Tv{z), 



(7-22) 



x{z,t) 




2w + 



7(1 — r) ST 7(1 + r) 



(7-23) 



1 — z z 1 + z 



u{z) = a + b - - 



~2w + 



27 



2-fz 



s 



(7-24) 
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The symmetry r — r is such that x{—z, — r) = x{z, r). 
We write: 

1 27 M 

Xo(z) = x(z,-l) = a + b 2w-\ , (7-25) 

2 1 — z z 

xi{z)=x{z,0) = a + b-l-2w + -^^, (7-26) 

2 I — z'^ 

1 u 2'y 

X2{z) = x{z, l)=a + b--~2w + - + (7-27) 

2 z 1 + z 

The coefficients w,s,7 are determined from the following system (we have written 
only the independent equations): 



3z suchthat X2{z) = Q and v{z) 
3z suchthat Xo{z) = b and v{z) 
3z suchthat xi{z) = a — ^ and v{z) 



2' 

1 

2' 



2' 



Those equations imply that: 

3w^ - w{2a + 2b-l) + b{a - 1) = 0, 

7 = 2w{b + 1 - w){a - w) = 2{w + l){b - w){a - w - 1), 
2s = -{2w + 1)(26 + 1 - 2w)(2a - 1 - 2w). 

And we have: 

3w = a + b h 27 + s. 

2 



(7-28) 



(7-29) 

(7-30) 
(7-31) 

(7-32) 



Envelope 
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Let s/27 = e. 
Parametrically 



2z^ 

One can check that it is the equation of the cardioid inscribed in the domain, see 
figiH 

The same domain can also be considered in the quantum case, but we don't do it 
here. 

7.3 The trapezoid 

Take = T, "D is the domain V{t^^^) = [0, 2T] at tmin = and '^'(tmax) = [f? ^] at 
^max = T. We choose the weights P = I and a ^ 1 with a constant in time. We study 
the classical case q = 1- 

Notice that 7^I'(tmin) = ^.N ^ N = #2^(tmax), therefore the initial matrix Mj^.^ 
is not fixed, it must have a potential Vt^^^ satisfying conditions eq. fl3-8p instead of 
eq. fl3-9p . but in fact it suffices to choose: 

= 0. (7-35) 

Let us look for an algebraic genus zero spectral curve of minimal degree. First, 
notice that we have no cusp condition at t^nm- We thus need to have only 2 poles, let 
us say at 2; = and z = 00, and the pole at 2; = disappears at t = T. Thus we write 
(with T = t/T): 

x{z, t) = 1 + c + rr + {ur + v)z . (7-36) 

z 

That means 

Xq{z) = x{z,{^) = 1 + c + vz , xi{z) = x{z,l) = l + c + r + {u + v)z. (7-37) 

Z 

and 

y,{z) = y,{z) = a Vj" ^^""^ (7-38) 
^ + {r + UZ + ^) 

The coefficients c, r, u, v can be determined by requiring that: 

• 3z such that xo{z) = and xi{z) = |. 

• 3z such that Xo{z) = 2 and Xi{z) = |. 

• In {yo{z)) ~ l/xo{z) at z ^ 00. 
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Figure 20: The envelope of the trapezoidal domain. The plots are for a — 
1,0.5,0.25,0.01,0.0001,200. 



The last condition implies m = and t = \ Then, the 2 other conditions imply 
c = and = (r — |)(r + |). 

Finally we find the following spectral curve: 

The envelope of the liquid region is: 

= 1 + 1 ^ ± ;;nrr vTi^ (7-42) 

z a + 1 a + 1 

this is a parabola shifted by a straight line, and tangent to at least two of the trapezoid 
boundaries. 

Notice that when a = 1, particles have the same probability to go upward hi — >■ 
hi + \ oi downward hi hi — ^, and the spectral curve is symmetric with respect to 
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X —>■ 2 — X. When a is very small, the probability to go downward is very small, and 
therefore almost all the particles are in the solid region going upward, and the liquid 
region becomes a narrow region around the line x = 1 — |. On the countrary, when a 
is very large, the probability to go downward is very large, and almost all the particles 
are in the solid region going downward, and the liquid region becomes a narrow region 
around the line x = 1 + |. See figj2Tl 




a=0.25 0=0.01 



Figure 21: Typical self-avoiding particles model, or typical plane partitions in the 
trapezoid. The tiling outside the liquid region is regular, it is frozen. For small a, the 
probability to go downward is very small, and therefore almost all the particles are in 
the solid region going upward, and the liquid region becomes a narrow region around 
the line x = 1 — |. 



7.4 The Plancherel law 

Let us choose a partition /i, and > n{^). we write: 

hiifi) = fii - t + N (7-43) 

we have: 

hiifi) > h2{n) > ...> hNifi) > (7-44) 

Consider the domain V, comprised between tmin = 0, and tmax = T, and such that at 
t = T the particles are at hi{T) = ^ + N — i (in some sense the boundary is the trivial 
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partition shifted by and at t = 0, the particles are at hi{0) = hi{ix) (the boundary 
is the partition fi). See fig 




Figure 22: The domain for the Plancherel law. At time t = t^ax = 
hi{traa.x) = ^ + N — i, and at time t = = tmin we have /li(tmm) = hi{fi). 



T, we have 



Let us call VrifJ') the plane partitions generating function in this domain: 

TT ,dlT\o=fJ.,dlT\T = 9 



Let us define the "Plancherel law" as the the limit T —>■ oo: 

Plancherel (/i) = V{fi) = "Pool/^)- 
It has been well known from a really long time [68j, that: 



P(^) 



and if g = 1, that reduces to the classical Plancherel law: 



nf=i/^.(/^)! 



(7-46) 



(7-47) 



(7-4J 



As a check of our matrix model approach, let us recover this classical result from the 
matrix model. 



63 



7.4.1 The classical Plancherel law q = 1 

As presented in section [3l the domain V is characterized by: 

• the matrix M^: 

Mt = I Id + diag(0, 1, 2, . . . , AT - 1), (7-49) 

• a potential at t = 0, which satisfies the conditions of eq. fl3-9p . we choose: 



e 



A'' 



Notice that we have: 



and thus: 



and 



U(y) = -T In (2 cosh |). (7-53) 



Theorem I3.H or more precisely its reduced version theorem 14.21 gives the relation- 
ship between the self-avoiding particles model partition function Pt(/^) and the matrix 
model: 

Z = |(iMoc^i?ie-^^«(^°)e-^^(^)e^^^(^--^°) 

i.e.: 

= C ^ 2 7-55 

And we can already guess, that in the large T limit, the role of the /ij(/i)'s in the 
spectral curve is going to be subleading, i.e., to large T leading order Z is going to be 
independent of hi{fi). 

In principle, our matrix model could be used to find the asymptotics of the 
Plancherel measure [601. 
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8 Obliged places and TSSCPPs 



So far, we have considered a self-avoiding particles model with defects, i.e. forbid- 
den places for the particles. Defects were introduced by choosing e~^*{x) = at the 
corresponding place {x,t). 

One could also be interested in constrained self-avoiding particles model, where we 
want to oblige some places to be visitted at some given times. This cannot be achieved 
directly by tuning the potentials Vf, but this can be achieved as follows. 

8.1 Obliged places 

We choose a potential Vt such that e~^'*^^) = 1 — 77 at the place {x,t), i.e. we enforce a 
probability 1—1] that the place {x, t) can be visited. In other words, the contribution 
of self-avoiding particles configurations such that one particle visits (x, t) will be pro- 
portional to 1 — r], and the contribution of self-avoiding particles configurations such 
that no particle visits (x, t) will be independent of 7] (notice that no more than one 
particle can visit {x,t)). Therefore the partition function Z is made of two terms: 

Z = Zo + (1 - r))Z^ (8-1) 

where Zi is the partition function of self-avoiding particles configurations which visit 
(x, t). We have: 



Z^^ - —Z 

drj 



2 i^T^Vtm) (8-2) 



r;=0 \ ' / ri=0 

In other words, Zi can be realized as some expectation value in our matrix integral. 
We write: 

-3) 



arj 



ri=0 



We have that: 

^^<TrA(Mt)>. (8-4) 

More generally, if we want to have several obliged places {xi,ti), i = 1, . . . , s, we 
introduce a function A{x, ti) for each of them. 



<Y[ TiA{Mt^,U)> . (8-5) 



Again, we have some arbitrariness in the choice of A{x,t). The only requirements 
are: 
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• and for t integer, tmin <t< tmax, we choose A{x,t) such that: 

^ ^ ^ ^ ^ ^ f 1 if 3i = {xi,ti) ^ ^ 

' 1 otherwise, m'D{t)/V, 

and arbitrary values everywhere else. 
A possible choice could be: 



sin {7t{x — Xi)) 

but many other choices could also be made. 



A{x,t) = y5t,u J 8-7 

■^-^ n [x — Xi) 



8.2 TSSCPP 

An application of what precedes is the partition function for counting TSSCPP (Totally 
symmetric self- complementary plane partitions) see fig [231 Counting TSSCPP has 
become a famous combinatorics problem, due to its link with alternating sign matrices, 
Razumov-Stroganoff conjecture, Hecke algebras and qKZ relations [211 H]. 

A TSSCPP configuration is completely determined by a partition of l/12th of the 
hexagon. 

In terms of a self-avoiding particles process, see fig. [2H this can be viewed as N 
self avoiding particles hi(t) jumping by ±|, and such that particle i has to follow a 
straight line after time t > N — i: 

hi{t) = N + - if t>N-i. (8-8) 
2 

In other words, we have a self-avoiding particles process with some obliged positions. 
Also, we don't fix the positions of the particles at time t = (although we could easily 
do it in our matrix model) . 

Notice that if we fix hi{N — 1) = N — i + at time t = N — 1, and if we oblige 
hi{N — z) = I (A^ — i) only at time t = N — i, then the self-avoiding particles process 
necessarily evolves in a way such that hi{t) = N — i + ^ ii t > N — i. In other words, 
it is sufficient to oblige only 1 position at each time: we oblige the position (y,t) at 
time t. 



8.2.1 The matrix model 

We choose for t = 0, . . . , A^ — 1: 

„ 2N-2-t 

A{x)=^ n (^-o-^') (8-9) 
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Figure 23: A TSSCPP, is a plane partition with all symmetries of the hexagon+ self 
complementarity (i.e. imagine that it is a pile of cubes within a big cube, then it must 
be equal to its complement). A TSSCPP configuration is completely determined by a 
partition of l/12th of the hexagon (the white region). 



where the constant Ct = l/t\{2N - 2 - 2t)\ is such that A(f ) = 1- 
We may also choose: 



For our matrix model, we choose the potential 14 = for all t (there is no defect). 
The TSSCPP partition function is thus realized by the matrix model: 




(8-10) 



where the constant Ct = ( 



1)7^! is such that A(f ) 




N-2 /N-1 



t=0 fi-l 



l) 



(8-11) 
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Figure 24: A plane partition of l/12th of the hexagon, can be reahzed as a self-avoiding 
particles process such that hi{t) = N- i + ^ift>N-i. 



where we integrate Mq, . . . , Mn-2 an i?i , . . . , -RAr_|, and Mjv-i = diag(^^^2~'"! ■ ■ i + 
— 1) is not integrated upon. All matrices are of size N x N. 

8.2.2 Examples 
• A^ = 2: 

2tsscpp(2) = I dModR^ Tr Ao(Mo) e"*'^''*^''^^ ^tri^, (m,-Mo) 

Mi = diag(^,^) , Ao(x) = ^(x-l)(x-2) (8-13) 



A^ = 3: 



Tr A/o ^Tr A/i 



^Tsscpp(3) = j dModMidRidRi Ti Aq{Mq) Ti Ai{Mi) q'' ""'^ q 

-trUiiRi) tri?i(Mi-A/o) -trUaiRs) tr R3 (Ma-Afi) /oi^\ 

e 1 ^ e ^ e 7 ^ e ^ (8-14) 

M2 = diag(l, 2, 3) , Ao(x) = ^(x- l)(x-2)(x-3)(x-4) (8-15) 

Ai(x) = ^(x-^)(x-^)(x-^) (8-16) 
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8.2.3 The spectral curve 

Since we have chosen all Vt = 0, we have the same spectral curve as for the Trapezoid 
of section I7.3[ Our partition function is computed as an expectation value. 



Wpp(iV) = V(lng)2^-2+^ Res... Res 

' ' OOl oojv 

g 

wi,'\zu...,ZN) A,{X,izi))...AMiXM{zN)) (8-17) 



8.2.4 Fixed position at time t=0 

If we want to fix the positions of particles at time t = 0, it suffices to introduce a 
potential 



e 



n2N~2f 



It would be interesting to relate those partition functions with fixed initial con- 
ditions for the particles, and see if we have some qKZ equations, and if they can be 
related to Alternating sign matrices partition functions. 

It was already noted [72J that there is a matrix model formulation for the 6-vertex 
model counting Alternating sign matrices. It would be interesting to compare the 
spectral curve of the 6-vertex matrix model, and the spectral curve of the matrix 
model we introduced here. 



9 Application to topological strings 

Regarding the applications to algebraic geometry and topological strings, our method 
allows to see that Gromov-Witten invariants of the toric Calabi-Yau 3-folds X = 
with branes (called the topological vertex), are the symplectic invariants of their mir- 
ror's spectral curve S^. In other words, we confirm that the "remodelling the B- model" 
proposal of BKMP [TT] holds for those toric CY 3-folds: 

BKMP's claim: GWg{X) = Fg{S^). (9-1) 

This claim is in the spirit of Dijkgraaf-Vafa conjecture, saying that B-model topological 
strings should be equivalent to some matrix model. Here, we have proved, that the 
topological vertex itself can be written as a matrix model. 

In order to prove BKMP for general toric CY 3 folds. It remains to find a "matrix- 
model way" to glue topological vertices together. 
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10 Conclusion 



We have introduced a new formulation of self-avoiding particles model in statistical 
physics, or lozenge tilings, dimer models, or plane partitions, in terms of matrix models. 
Then we briefly discussed how to apply the rich technology of matrix models. 

However, this is only a starting point, it is only a proposed framework, and we 
expect that exploiting the immense knowledge of matrix integrals developped since 
Wigner-Dyson-Mehta, we can find new consequences for self-avoiding particles model. 

It should be possible to look at many extensions of this approach, for instance taking 
appropriate limits, it should be possible to apply this formalism to Dyson processes, 
or vicious walkers, and recover many results in the matrix model framework. 

As an application to algebraic geometry, we have only found the "topological ver- 
tex", i.e. the building block for computing Gromov-Witten invariants of all toric CY 
3-folds, but unfortunately, it is at the present time, not known how to perform the 
gluing of vertices, in the language of symplectic invariants of [27] . This question needs 
to be addressed. 
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A Gamma-function and q-product 

A.l Gamma- function 

The Gamma-function T{x) is such that: 




(1-1) 
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If X = n + 1 is a positive integer: 

r(n+l) = n! =n(n-l)(n-2)... (1-2) 
r has poles at all negative integers: 

r{-n + e)r^^-^ (1-3) 



We have: 



r{x)r{-x) = . J . (1-4) 

X sm {ttx) 



Stirling formula: 



-0(0;) ~ Inx — r — (1-5) 



^ ' n=l 

where 5^ are the Bernouilli numbers. 

^1 = -^ , B2^l , ^4 = -— (1-6) 
2 ' 6 30 ^ ^ 

A. 2 q-product 

The q-product is defined as an infinite product: 

oo 

The function rg(a;) shares many similarities with the function T{x), by replacing inte- 
gers with g— numbers. It is a quantum F-function. g— numbers are: 



N ^q^ -q^ (1-8) 



They tend to usual numbers in the limit g — > 1. 
If n is an integer we have: 



n-l 

9{ql = , g'{ql=g{l)q-'^ llim]=g{l)q-'^ [n-l]\ (1-9) 



m=l 

Stirling formula at small In q: 



ln^(a;) = ^ E ^, " (1^?)" ^^2-n(l/a:) (1-10) 



n=0 



^ ^(x) ^ n! 
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where Lin is the Polylogarithme: 



OO u 

X 



:i-i2) 



Lii{x) = — In (1 — x) 

That is: 

xg'{x) 1 
g{x) Ing 



Lio(x) 



fc=i 

X 

1 — X 



1 



= -L^„_i(x) (1-13) 

X 



n=l 



B Introduction to symplectic invariants 

Here, we briefly recall the definition of symplectic invariants, but we refer the reader 
to [271 [2H] for more details. 

They were first introduced in [2H[T5] and further formalized in [27], as a solution of 
the so called "topological expansion" of matrix integrals. But then, in [27], they were 
defined as algebraic quantities for spectral curves, independently of the existence of an 
underlying matrix model. 

B.l Definition of symplectic invariants 

Consider a spectral curve S = (£, x, y), where £ is a compact Riemann surface, and x 
and y are two analytical functions, such that dx and dy are meromorphic forms on C 
The branchpoints are defined as the zeroes of dx: 

dx{ai) = (2-1) 

We assume that the spectral curve is regular, i.e. we have a finite number of branch- 
points, and they are simple branchpoints, which means that each Oj is a simple zero 
of dx, and dy{a.j) ^ 0. This is equivalent to say that near Oj, we have a square root 
behaviour: 

y{z) ~ y{ai) + C ^x{z)~x{a{). (2-2) 

Since we have a square root branchpoint near a^, this means that in the vicinity of 
Oj, there exists z ^ z on the other branch, i.e. x{z) = x{z): 

nearoi, 31^7^2;, x{z) = x{z). (2-3) 
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z is called the conjugate of z near a^. It is defined only locally near branchpoints, and 
it is not necessarily defined globally. 

On a Riemann surface £, there exists a "Bergman kernel", i.e. a symmetric mero- 
morphic 2-form, having a double pole, with no residue, on the diagonal, i.e. which 
behaves like: 

„ / \ d,z^ dzo ^ , s / . s 

5^1,^2 ~ , 2-4 

21^22 [Zl - 

where z is any local parameter. 

The Bergman kernel is unique if we also fix its periods §^ B, but this is not nec- 
essary for the definition of symplectic invariants. Let us assume that we have choosen 
one Bergman kernel B{zi, Z2). 

Definition B.l We define the following n— forms: 

Wi:°\z) = ~y{z)dx{z) (2-5) 

lyf (^1,^2) = 5(^1,^2) (2-6) 

and by recursion for 2 — 2g — n < and n > 1, and if we denote collectively J = 
{z2, . . . , Zn}-' 

Wjf\z,,J) = VResif(^i,^) [Wif/)(^,^,J) 

i 

+ E E ^) ^n-wi-^^ JM)]- (2-7) 

h=0 ICJ 

where: 

For n = 0, and g >2, we define: 

^ i 

where $ is any function such that d^ = ydx. 

The Wi^^'s defined in this way are always symmetric meromorphic n-forms, having 
poles only at branchpoints (except Wi^^ and W^2°^). And Fg does not depend on the 
choice of integration constant for $. 

There is also a definition for Fq and Fi, but we refer the reader to p7j. 
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B.2 Symplectic invariants of genus zero spectral curves 

An interesting example is when L is the complex plane C. x and y are functions of a 
complex variable z. In that case we have: 

For example we have: 

wt\z,,Z2,z,) = J2^„^Jy,^^^^ ^,^_a,Y%^X)hz,-a.?- ^^"^^^ 

B.3 Examples of spectral curves and their symplectic invari- 
ants 

Let us give a few examples of spectral curves: 

• The spectral curve Swp defined by the functions x{z) = z'^,y{z) = ^sin(27r2;), 
i.e. y = 2^ sin (27ry^) appears in the Weil-Petersson volumes of moduli spaces, its Fg's 
are the Weil-Petersson volumes: Fg{Swp) = Vol(A^c,), see [29]. 

• The spectral curve ^Airy defined by y^—x = 0, i.e. y = y/x is associated to the Airy 
kernel law, and to the universal behavior of extreme eigenvalues distribution, i.e. to 
the Tracy- Widom law [70]. The Fg's of the Airy curve are equal to zero: Fg(5Airy) = 0. 

• The spectral curve y = ^yPo\2m+iix), where Pol2m+i is a polynomial of degree 
2m + 1, is associated to the (2m + 1, 2) minimal model in conformal field theory (with 
central charge c = 1 — 3 ^2^1+1 )' ^ reduction of the KdV hierarchy. Its F^'s 
are related to the KdV Tau-function, and they can be computed by the asymptotic 
expansion of the solution of a Painleve type equation (in fact the m + 1*^* Gelfand-Dikii 
equation Rm+i{u{t)) = t). The case m = i.e. the (1,2) minimal model is the Airy 
case. The case m = 1 i.e. the (3, 2) minimal model is called pure gravity, and the case 
m = 2 i.e. the (5, 2) minimal model is called Lee- Yang singularity. 

• Spectral curves of type Pol(e^,e^) = appear in topological strings, where Pol is 
a polynomial. More precisely, consider a Toric Calabi-Yau 3-fold X. Through mirror 
symmetry [ID] , it has a mirror X, which is also a toric Calabi-Yau 3- fold, and which is 
given by a submanifold of of equation H{e^, e^) = uv where if is a polynomial. The 
spectral curve S^, is defined as the singular locus of X, which satisfies the equation 
H{e^,ey) = 0. The "remodelling the B-model" conjecture of BKMP [11], is that the 
generating function of Gromov-Witten invariants GWg{X) of genus g of X, is equal to 
the symplectic invariant Fg{Sj^) of the spectral curve of the mirror: 

GWg{X) = Fg{S^) (2-12) 
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This conjecture was proved in few cases, and in particular for the family of Hirzebruch 
manifolds Xp = 0{—p) © 0(p — 2) ^ (which include the connifold), see [26]. It 
was also partially proved for SU{n) Seiberg-Witten theories |16], and for 2* geometries 

B.4 Some properties of symplectic invariants 

Let us now give a few properties of symplectic invariants. We refer the reader to [271 [21] 
for the general theory, and for proofs and details. 

Properties: 

• Homogeneity: If we rescale y —>■ Xy, which we note XS = (x, Xy), we have: 

Fg{XS) = X^-^^ Fg{S) , W^'^\XS) = X^~^^-''Wjf\S). (2-13) 

In particular that implies Fg{—S) = Fg{S). 

• Symplectic invariance: if two spectral curves S = {y{x)} and S = {y{x)} are 
symplectically equivalent, i.e. if there is a map from C x C to C x C which sends one 
spectral curve to the other and conserves the symplectic form dx Ady = dx A dy, then: 

FgiS) = FgiS) (2-14) 

For instance we can change y y + R{x) where R{x) is a rational function, or [y — >■ 
x, X — > —y), or {x — > \iax,y — > xy), ... etc, without changing Fg. 

The W!f^ with n > 1 are not symplectic invariants, they change under symplectic 
transformations, but they change in a covariant way. For example wi^\S) - wi'^\S) 
is an exact form. 

• Variations: Consider an infinitesimal deformation of the spectral curve (x, y) — > 
{x + e6x,y + e6y), such that Q = 6y dx — 6xdy is a meromorphic differential form on 
the spectral curve. Any meromorphic form Q is dual to a cycle Q* on the curve, the 
duality pairing is realized through the Bergman kernel: 

n{z) = (f B{z,z') (2-15) 
Jn* 

Then, the infinitesimal variation of Fg and W^f '' is given by: 

^ Wjf\z,, ...,z„)= £ Wi%iz,, . . . , z„, z') (2-16) 

The infinitesimal variation of Fg is the case n = 0: ^ = ^j^, Wi^\z'). 

• Limits: Consider a one parameter family of spectral curves S{t), such that at 
t = tc, the spectral curve S(tc) becomes singular. Consider its blow up: S(t) ~ 
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{t — tcY Sc + o((t — tc)*^); where the exponent v is chosen such that the curve Sc is 
regular. Then we have: 

Fg{S{t)) ^{t- t,)(2-23). ^^(5^) + _ (2-17) 
and more generally 

This theorem is very useful, it gives the asymptotic behaviors. For instance, if we 
zoom near a regular branch point, the spectral curve always behaves like y = \fx^ 
i.e. Sc = SAiryi and we find the Airy law, and Tracy- Widom. Near an algebraic cusp 
singularity y ~ x^/'', we find the spectral curve Sc = S^p^q) of the (p, q) minimal model 
of conformal field theories (of central charge c = 1 — 6 ), which is a reduction of 
the KP hierarchy. 

Thes are just a few of the properties satisfied by the F^'s. See |27j for more. For 
instance there are also modular properties and holomorphic anomaly equations, Hirota 
equations, and diagrammatic representations. 

C Solution of matrix models 

Let us present a brief review about matrix models. 

C.l Generalities on the solution of matrix models 

Consider a matrix integral of type: 

Z= f f\dM, ^-Q^[T.Uv.m)+j:.c,M,M,+,] (3_^) 

where HnICi) is the set of normal matrices having their eigenvalues on contour Cj, and 
where Mp+i is not integrated upon. 

If one wants to find a large Q expansion of In Z, the answer, is that one has first 
to compute the "spectral curve" S (defined below), then, in p3], it is proved that: 

Theorem C.l If S is the spectral curve associated to potentials Vi and paths Ci, and 
if Z has a large Q expansion of the form 

oo 

\TiZ = Y,Q^~^'Fg{S) (3-2) 

9=0 

then, the coefficient Fg = Fg{S) is the symplectic invariant of degree 2 — 2g of the 
spectral curve S. Symplectic invariants Fg{S) of a spectral curve S were defined and 
introduced in l2l^ . 
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We shall explain the meaning of this theorem below in further details. We shall 
explain how to compute the spectral curve 5 of a matrix model, and how to compute its 
symplectic invariants Fg. We recall the definition of symplectic invariants in appendix 
11 

Let us just mention that finding the spectral curve is rather automatic, it is mostly 
an algebraic task, and for most examples the spectral curve is a rather simple object. 
Here, we have the the spectral curves of Kenyon-Okounkov-Sheffield |48] . 

Then, computing the symplectic invariants of a spectral curve, is rather easy. The 
definition of symplectic invariants [27] involves computing residues of rather standard 
functions (exponentials, logs, rational functions, ...), and can be completely autom- 
atized, see appendix [Bl It is really an efficient method. For example Fi can often 
be computed by hand. Moreover, symplectic invariants satisfy many properties, which 
make them really convenient to use, for example there are formulae for computing their 
derivatives with respect to any parameter, and formulae for finding their limits near 
singularities. 

C.2 Generalities about loop equations 

The theorem IC.ll for the chain of matrices was proved from loop equations [33] . 

Loop equations are obtained by integration by parts in the matrix integral, or 
equivalently, by writing that an integral is invariant under change of variable. For 
example, assume that we make an infinitesimal change of variable in eq. ( 13-11) : 

Ms ^ M3 + eM|M| + 0{t^) (3-3) 

one finds the Jacobian to first order in e: 

dM^ dM^ [l + e( Tr M| Tr M| + Tr M^Mg Tr M3 + Tr M|M| Tr Id) + 0{e^)^ (3-4) 

Writing that the integral is invariant to first order in e implies: 

< Tr M| Tr M| + Tr M^Mg Tr M3 + Tr M|M| Tr Id > 
= Q < Tr Vl{M^)MlMl + C2 Tr M|M| + C3 Tr M^M^Mi > (3-5) 

Similarly, by considering other appropriate changes of variables, one can write a fam- 
ily of relationships of the type eq. ( I3-5P between expectation values, this was done 
systematically in [33J. 

Instead of considering expectation values of powers < Tr >, is is more con- 
venient to group them in a formal generating function Wi{x) =< Tr ^^^-^y >= 
^'^fc+i''^ - This leads to introduce: 

.... x„) =< TV ^ ^ ... TV ^ (3-6) 
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(where the subscript < . >c means the cumulant, or connected part). Since we made 
the hypothesis that there exists a large Q expansion we write: 

oo 

W^{x^, . . . , = 5^ Q'-'^-" Wi^Kxi, . . . , x„). (3-7) 

9=0 

The problem consists in finding a family of relationships among those quantities, 
and which allow to find the solution, i.e. compute all of them. This was done in [531 E2]- 

The result of [33], is that the Wn\xi, . . . , Xn) = Wn\xi, . . . , x„) dxi . . . dx^s are 
the n-forms of theorem \C.1\ i.e. the Wn^^^s of [27j, and where the spectral curve 
S = (x, y) is y = —W^^\x). 

In some sense, (x) can be viewed as the " large Q" limit of the resolvent of the 
matrix Mi (the first matrix in the chain eq. fl3-ip ). it is often called the "equilibrium 
distribution of eigenvalues of Mi", but one should take this denomination with some 
care. Indeed, W^f ^ can be shown to be the "weak" large Q limit of the resolvent, 
for some potentials (in particular in the potentials don't depend on Q), and some 
integration paths Cj, but this is probably wrong in general. In fact, in our case, the 
potentials do depend on Q, and saying that the spectral curve is the "large Q" limit 
of the resolvant is slightly wrong (and at least one needs to make precise the notion of 
large Q limit). 

The only correct definition of the spectral curve, which corresponds to theorem lC.il 
IS that y = Wi^^^ (x), is the first term in the formal large Q expansion, doing as if the 
potentials were independent of Q. 

Working to order g = in the expansion eq. fl3-7p within the loop equations, is 
formally equivalent to replacing expectation values of product of traces, by product 
of expectation values of each trace (indeed, the connected part comes with a factor 
Q2-2g-n^ whcrcas the non-connected term, i.e. the factorized term comes with a factor 
Q2-2g+n^ i.e. all terms which are not disconnected don't contribute to the highest power 
of Q): 

< Tt AiTr A2 . . . Tt Ak > — ^ < Tr Ai >< Tr > . . . < Tr > (3-8) 

This formal manipulation allows to find an equation which determines the spectral 
curve Wi^\x). 

C.3 Spectral curve of the chain of matrices 

Here, we shall just give a "ready to use" recipe of how to find the spectral curve for 
a chain of matrices matrix model of type eq. (13-11) . This recipe is extracted from [33] . 
and it is really technical to explain how to obtain it. However, it is easy to use. 
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C.4 General case 

The spectral curve of a matrix integral of type eq. (13-ip : 

Z = [ rTrfM,e-«^E'-^«(*^^)+^-^»^«*^»+^l (3-9) 

is characterized by a set (see [321 [33]) of p + 2 analytical functions of a variable z {z 
belongs to a Riemann surface C). There is one such analytical function for each matrix 
Mi, i = 1, . . . ,p + 1 of the chain, plus one additional function at the end of the chain. 
Let us call them: 

X,{z) , t = 0,...,p + l. (3-10) 
Those functions are completely determined by the following system of equations: 

\fi = 2,...,p, q_iX,_i(^) + qX,+i(z) + V;'(X,(^)) =0 (3-11) 

and 

wl''\M^)) = M^) = viiM^)) - ciM^), (3-12) 

together with the conditions: 

• Xp{z) has simple poles at the values of z such that Xp+i{z) is an eigenvalue of 
Mp_|_i. Let us call d, the value of z such that Xp_|_i(2;) = Aj the i^^ eigenvalue of Mp+i. 
Then we have in the vicinity of z ^ Q: 

Xp{z) ~ \ (3-13) 

-^p+i{z) — Ai 

• The function W^^\x), is analytical outside of the cuts [0^,64], the endpoints of 
the cuts being zeroes of dXi. Near Xi ^ 00 (only in the physical sheet), we have: 

Mz)-^ + 0{l/Uzr) (3-14) 

• the genus of the Riemann surface C, and the period integrals £^ XodXi on non- 
contractible cycles Ai of C, are related to the choice of contours Ci in the matrix integral 
eq. f l3-9p . The following quantity is called "filling fraction" 



1 

2i7T 



XodXi, (3-15) 



A 



and it is such that Qei is the number of eigenvalues of Mi contained in the domain 
surrounded by the projection in C of the cycle Ai by the application Xi. 

More generally, ^ §^ Xj^idXj is the number of eigenvalues of Mj contained in 
the projection in C of the cycle Ai by the application Xj. 
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Finding the spectral curve for arbitrary potentials Vi and arbitrary contours Ci 
can be extremely tedious. However, here we are interested in formal matrix integrals, 
which are defined as formal power series expansions in some formal parameter, and in 
particular, the spectral curve we are looking for, must also be found as a formal power 
series. In particular, the genus is the genus of the spectral curve when we send the 
formal parameter to 0, and in general that simplifies the computation of the spectral 
curve a lot. Very often the genus is in fact zero. 

C.4.1 Topological expansion 

When we have determined the spectral curve S, i.e. the two functions x{z) = Xi{z) 
and y{z) = Xo{z): 

S = ix,y) = Si,o = (Xi,Xo) , xiz) = X,{z), y{z) = Xo(^), (3-16) 
we have: 

oo 

lnZ = 5^F,(5). (3-17) 

9=0 

It is then useful to use the symplectic invariance of the Fg^s. 
Consider the following spectral curves: 

Si,j = iXi,Xj) (3-18) 

Due to eq. (13-1 ip . we have: 

CidXi+i A dXi = -Ci-idXi^i A dXi = Ci^idXi A dXi_i (3-19) 

i.e. all the spectral curves are symplectically equivalent: 

CiSi^i^i = Cj_nSj_i^j = — Cj_nSj^j_i (3-20) 

and therefore, the theorem of symplectic invariance of the F^'s [271 [30] gives Vi: 

c^-'' F,(5,,+i) = c^-^^ F,(5,_m) = cl',' F,(5,,_0 (3-21) 

In other words, the F^'s can be computed with the spectral curve S = {x, y) where x 
and y are any two consecutive we don't need to choose the pair Xi,Xo. 

C.4.2 Densities and correlation functions 

For the general chain of matrices, the wlf^^s were defined as the Stieltjes transforms of 
the density correlation functions of the first matrix Mi in the chain, i.e. the resolvents 
eq. fl3-6p . The densities can be recovered by taking the discontinuities: 

p{h) = J2Q'-'' P^'\h) = (^J2^{h-h,)^ (3-22) 
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where 

(h) = (wi'^ {h - lO) - Wi'^ {h + zO) V (3-23) 

We recall that the spectral curve used to compute w!f^ here, is the spectral curve 
Sifi = {Xi,Xq), and the wjf^'^s are not invariant under symplectic transformations. 

However, just by applying eq. fl2-16p . one can see that, for any matrix Mj of the 
chain, one has that the density of eigenvalues of Mj: 

p,{h)=Y,Q'-''pf{h) = TrSih-M,)^ (3-24) 

is given by: 

pf{h) = \ fwf'\h - ^0) - Wi'\h + zO)) . (3-25) 

Therefore, the density of eienvalues of each matrix of the chain, can be computed 
to all orders in the large Q expansion, using the W^^\ defined 

In the general chain of matrices, the support of the densities must be related to the 
integration paths Cj's: One defines the support of densities as the loci where densities 
are real. The supports of densities, together with their analytical continuations, must 
be homologicaly equivalent to the integration paths Cj's. And saying that the paths 
Ci are steepest descent paths, is more or less equivalent to saying that the densities pi 
must be real and positive on their supports, and i times the analytical continuation 
of the densities, along the analytical continuations of the supports, must be real and 
strictly positive (stability condition). This is a highly non trivial condition in general. 

Fortunately, for formal matrix integrals, defined as formal power series in a formal 
parameter, the small limit of the parameter often gives a very simple matrix model, 
with a very simple spectral curve, and it is often easy to check positivity to leading 
order, and then, that the formal power series expansion does not destroy the supports 
and positivity. 
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